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' Abstract. We prove existence and uniqueness of stochastic representations for solutions to elliptic 

and parabolic boundary value and obstacle problems associated with a degenerate Markov diffusion 
. process. In particular, our article focuses on the Heston stochastic volatility process, which is widely 

' used as an asset price model in mathematical finance and a paradigm for a degenerate diffusion 

process where the degeneracy in the diffusion coefficient is proportional to the square root of the 
distance to the boundary of the half-plane. The generator of this process with killing, called the 
elliptic Heston operator, is a second-order, degenerate, elliptic partial differential operator whose 
Qi^ , coefficients have linear growth in the spatial variables and where the degeneracy in the operator 

(-H ' symbol is proportional to the distance to the boundary of the half-plane. In mathematical finance, 

solutions to terminal/boundary value or obstacle problems for the parabolic Heston operator cor- 
respond to value functions for American- style options on the underlying asset. 
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1. Introduction 

Since its discovery by Mark Kac [29], inspired in turn by the doctoral dissertation of Richard 
Feynman |19j . the Feynman-Kac (or stochastic representation) formula has provided a link between 
probability theory and partial differential equations which has steadily deepened and developed 
during the intervening years. Moreover, judging by continuing interest in its applications to mathe- 
matical finance |31j and mathematical physics [37^ |4T] , including non-linear parabolic equations [7] , 
this trend shows no sign of abating. However, while stochastic representation formulae for solutions 
to linear, second-order elliptic and parabolic boundary and obstacle problems are well established 
when the generator, —A, of the Markov stochastic process is strictly elliptic [H [Ml [Ml |38] in the 
sense of |25l p. 31], the literature is far less complete when A is degenerate elliptic, that is, only has 
a non-negative definite characteristic form in the sense of [39], and its coefficients are unbounded. 

In this article, we prove stochastic representation formulae for solutions to an elliptic boundary 
value problem, 

Au = f ouG, (1.1) 

and an elliptic obstacle problem, 

min{An - /, n - ^} = on 0, (1.2) 

respectively, subject to a partial Dirichlet boundary condition, 

u = g on Fi. (1-3) 

Here, the subset ^ EI is a (possibly unbounded) domain (connected, open subset) in the open 
upper half-space EI := M.'^~^ x (0, oo) (where d > 2), Fi = dff R EI is the portion of the boundary, 
d^, of which lies in EI, / : ^ — )• M is a source function, the function g : Fi — )• M prescribes 
a Dirichlet boundary condition along Fi and ^/^ : ^ U Fi — )• M is an obstacle function which is 
compatible with g in the sense that 

V'<5 onFi, (1.4) 
while A is an elliptic differential operator on which is degenerate along the interior, Fq, of 
dM n 80' and may have unbounded coefficients. We require Fq to be non-empty throughout this 
article as, otherwise, if is bounded (and the coefficients of A are, say, continuous on 0'), then 
standard results apply [5l [2^ [30l [38] . However, an additional boundary condition is not necessarily 
prescribed along Fq. Rather, we shall see that our stochastic representation formulae will provide 
the unique solutions to (jl.ip or (jl.2p , together with ()1.3p , when we seek solutions which are suitably 
smooth up to the boundary portion Fq, a property which is guaranteed when the solutions lie in 
certain weighted Holder spaces (by analogy with [llj). or replace the boundary condition (jl.3p 
with the full Dirichlet condition, 

u = g on 80, (1.5) 
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in which case the solutions are not guaranteed to be any more than continuous up to Fq and 
ip : is now required to be compatible with g in the sense that, 

t/j <g on diff. (1.6) 

We also prove stochastic representation formulae for solutions to a parabolic terminal/boundary 
value problem, 

-ut + Au = f onQ, (1.7) 

and a parabolic obstacle problem, 

mm{—ut + Au — f,u — 'ip} = onQ, (1-8) 
respectively, subject to the partial terminal/boundary condition, 

u = g on g^Q. (1.9) 
Here, we define Q := (0, T) x 0', where < T < oo, and define 

g^Q := (0,T) X riU{r} X (^UTi), (1.10) 

to be a subset of the parabolic boundary of Q, and now assume given a source function / : Q — t- M, 
a Dirichlet boundary data function g : S^Q — t- M, and an obstacle function ip : QU Q^Q — > M which 
is compatible with g in the sense that, 

'il'<g ong^Q. (1.11) 

Just as in the elliptic case, we shall either consider solutions which are suitably smooth up to 
(0,T) X Tq, but impose no explicit Dirichlet boundary condition along (0, T) x Tq, or replace the 
boundary condition in (jl.9p with the full Dirichlet condition 

u = g ondQ, (1.12) 

where 

gg := (o,r) X 9^u{r} X ^, (1.13) 

is the full parabolic boundary of Q, in which case the solutions are not guaranteed to be any more 
than continuous up to (0, T) x Tq and ip : QU dQ — )• M is now compatible with g in the sense that 

ij<g on gg. (1.14) 

Before giving a detailed account of our main results, we summarize a few applications. 



1.1. Applications. In mathematical finance, a solution, u, to the elliptic obstacle problem (jl.2p . 
(jl.3p . when / = 0, can be interpreted as the value function for a perpetual American- style option 
with payoff function given by the obstacle function, ip, while a solution, u, to the corresponding 
parabolic obstacle problem (II. Sp . ()1.9|) . when / = 0, can be interpreted as the value function 
for a finite-maturity American-style option with payoff function given by a terminal condition 
function, h = g{T, •) : ^ — t- M, which typically coincides on {T} x 0' with the obstacle function, 
■ip. For example, in the case of an American-style put option, one chooses il:{x,y) = {E — e^)"*", 
\/{x,y) G ^, where > is a positive constant. While solutions to (II. ip . (II. Sp do not have an 
immediate interpretation in mathematical finance, a solution, u, to the corresponding parabolic 
terminal/boundary value problem (II. 7|) . (II. 9p . when / = 0, can be interpreted as the value function 
for a European- style option with payoff function given by the terminal condition function, h. For 
example, in the case of a European-style put option, one chooses h{x,y) = {E — e^)"*", V(x,y) G ^. 

Stochastic representation formulae underly Monte Carlo methods of numerical computation of 
value functions for option pricing in mathematical finance [26]. As is well-known to practitioners. 
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the question of Monte Carlo simulation of solutions to the Heston stochastic differential equation 
is especially delicate [3l[36]. We hope that our article sheds further light on these issues. 

1.2. Summary of main results. In this article, we set d = 2 and choose —A to be the generator 
of the two-dimensional Heston stochastic volatility process with killing rate r |27] . a degenerate 
diffusion process well known in mathematical finance, 

- At; := I {v^^ + 2pav^y + a'^Vyy) + {r - q - y/2)v^ + k(0 - y)vy -rv, v£ C~(M). (1.15) 

Nonetheless, we expect that many of our results would extend to a much broader class of degenerate 
Markov processes and we shall address such questions elsewhere. Throughout this article, the 
coefficients of A are required to obey 

Assumption 1.1 (EUipticity condition for the Heston operator coefficients). The coefficients defin- 
ing A in (jl.lSp are constants obeying 

C7/0,-l</3<l, (1.16) 

and K > 0, 6* > 0, ancfl q', r G M. 

Let (fl,^,F, Q) be a filtered probability space satisfying the usual conditions, where F = 
{.^(s)}s>o is the Q-completion of the natural filtration of {W{s))s>q-, and iW{s))s>o is a stan- 
dard Brownian motion with values in M^. For < t < T < oo, let ^^t,T denote the set of F-stopping 
times with values in [t, T]. Let (X*'^'^(s), y*'^(s))s>( denote a continuous version of the strong 
solution to the Heston stochastic differential equation 

dX{s)={r-q-^^^ ds + y/Y{s)dWi{s), s > t, 

dY{s) = K(i9-y(s)) ds + ay^Y{^ (^pdWiis) + y^l - p"^ dWiis)^ , s > t, ^^'^^"^ 
iX{t),Y{t)) = ix,y), 

which exists by Corollary 12.81 where the coefficients are as in Assumption 11.11 For brevity, we 
sometimes denote z = {x,y) and {Z^'^{s))s>t = {X^'^'^ {s),Y^'y {s))s>t- We omit the superscripts 
(t, z) and (t, X, y) when the initial condition is clear from the context, or we omit the superscript t 
when t = 0. We let 

/5:=^, (1.18) 

M:=^, (1.19) 
denote the Feller parameters associated with the Heston process. 

1.2.1. Existence and uniqueness of solutions to elliptic boundary value problems. For an integer 
k > 0, we let C^[iff) denote the vector space of functions whose derivatives up to order k are 
continuous on ^ and let C^{d') denote the Banach space of functions whose derivatives up to order 
k are uniformly continuous and bounded on ^ [21 §1.25 & §1-26]. If T ^ is a relatively open 
set, we let C]^^(i^ur) denote the vector space of functions, u, such that, for any precompact open 
subset f/ (s ^ U T, we have u £ C^{tj). 
We shall often appeal to the following 



^We impose additional conditions, such as q > 0, r > 0, or r > 0, depending on the problem under consideration; 
we only require that g > when deriving the supermartingale property in Lemma l2.11l IT}, a property used only in 
the elliptic case. 
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Hypothesis 1.2 (Growth condition). If w is a function then, for all {x, y) in its domain of definition, 

\vix, y)\ < C{1 + e^^^^y + e^-^^x)^ (1^20) 
where C > 0, < Mi < min {r/ (k??) , fx}, and M2 G [0, 1). 
Let [/ ^ H be an open set. We denote 

:= inf {s>t: Z^'\s) ^ U} , (1.21) 

and we let 

inf {s>t: Z*''{s) (^Uu{Undm)}. (1.22) 

r^^ = i/^^. We also have that r^^ = z/^^ 
in this case the process Z^'^ does not reach the boundary dM, by Lemma I2.1UI ([T]). By |38[ p. 
117], both and z/^^ are stopping times with respect to the filtration F, since F is assumed to 
satisfy the usual conditions. When the initial condition, {t,z), is clear from the context, we omit 
the superscripts in the preceding definitions ()1.2ip and (ll.22|) of the stopping times. Also, when 
t = 0, we omit the superscript t in the preceding definitions. 

Theorem 1.3 (Uniqueness of solutions to the elliptic boundary value problem). Let r > 0, q > 0, 

and f belong to C{&') and obey the growth condition (jl.20p on G . Then 

(1) If 13 > I, assume g £ Cioc(ri) obeys HfLlU^ . Let 

nGCioc(^uri)nc2(^) 

be a solution to the elliptic boundary value problem (|l.ip . (II. 3p and which obeys (ll.20p on 
G . Then, u = u* on 0" L)Ti, where 



t.z 
Vj} -.= L 

Notice that if C7 R dM = 0, then rjf = ujf . We also have that rjf = vlf when /? > 1, because 



u*{z) :=E^ [e-'''"^<7(^(r^))l{r^<oo}] +n 



'f{Z{s))ds 



(1.23) 



where Tg is defined by (jl.2ip . for all z £ U Fi. 
(2) IfO<l3< I, assume g G Cioc(<9^) obeys (fT:20]l on dff, and let u £ Ciod^) nC^ {^) be 
a solution to the elliptic boundary value problem (II. ip . (jl.Sp and which obeys (ll.20p on G . 
Then, u = u* on G , where u* is given by (ll.23p . 

Following m Definition 2.2], we let C]'/^^(^UFo) denote the subspace of C^'^'iG) nCl^iGuFo) 
consisting of functions, u, such that, for any precompact open subset U (£ G U Fq, 

sup \u{x, y)\ + \Du{x,y)\ + \yD'^u{x,y)\ < Qo, (1-24) 

{x,y)eU 

where Du denotes the gradient and D^u the Hessian matrix of u, defined Lebesgue-a.e. on G. 

Remark 1.4 (On the definition of C^\qj,(^UFo)). The definition of the function space cI'1^^{G'UTq) 
can be relaxed in our article from that of [El Definition 2.2]. It is enough to say that a function u 
belongs to Clf^^{ff U Fq) if u G C^'^{G) n Cioc{G U Fq) and (fL2ip holds on any precompact open 
subset [/ (s ^'uFq. 

Theorem 1.5 (Uniqueness of solutions to the elliptic boundary value problem (jl.ip . (]1.3p . when 
< ;S < 1). Lei r > 0, g > 0, < /3 < 1, and let f be as in Theorem HOI Let g G Cioc(Fi) obey 
(ll.20p on Fi and suppose that 

u G c,U^ u Fi) n c\G) n c]i^{G u Fq) 
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is a solution to the elliptic boundary value problem (jl.ip . (I1.3p which obeys ()1.20p on ff. Then, 
u = u* on ffUTi, where u* is given by 



u*(z) :=E^ [e-^'''^g{Z{uff))l{,^^^y] + E| 
and Vff is defined by ()1.22p . for all z ^ ^ U Fi. 



r e-^'f{Z{s))ds 
Jo 



(1.25) 



Remark 1.6 (Existence and uniqueness of strong solutions in weighted Sobolev spaces to the ehiptic 
boundary value problem). Existence and uniqueness of strong solutions in weighted Sobolev spaces 
to problem (ll.ip with boundary condition (I1.3P along Fi, for all /3 > 0, is proved in [U Theorem 

I. 18], and Holder continuity of such solutions up to Fq is proved in |151 Theorem 1.10]. 

Remark 1.7 (Comparison of uniqueness results). To obtain uniqueness of solutions to the elliptic 
boundary value problem (jl.ip with boundary condition p.3p only specified along Fi, we need to 
assume the stronger regularity hypothesis 

u G Cioc(^ u Fi) n c2(^) n c];L(^ u Fq) 

when < /3 < 1, while the regularity assumption 

nE C7ioc(^UFi)nC2(^) 

suffices when /? > 1. The analogous comments apply to the elliptic obstacle problems described in 
Theorems 11.121 and I1.13[ the parabolic terminal/boundary value problems described in Theorems 

II. 161 and I1.18[ and the parabolic obstacle value problems described in Theorems 11.271 and 11.281 

For a G (0,1), we let C'^"^"(€?') denote the subspace of C^{ff) consisting of functions whose 
derivatives up to order k are locally a- Holder continuous on (in the sense of \25\ p. 52]) and let 
Qk+a^^^ denote the subspace of C^{S') consisting of functions whose derivatives up to order k are 
uniformly a-Holder continuous on ^ |25l p. 52], [21 §1.27]. If T ^ dff is a relatively open set, we 
let Cj^|'*(^UT) denote the vector space of functions, n, such that, for any precompact open subset 
[/ ^ ^ U T, we have u G C'=+"(C7). 

We have the following result concerning existence of solutions to the elliptic boundary value 
problem with traditional regularity on 0'. 

Theorem 1.8 (Existence of solutions to the elliptic boundary value problem with continuous 
Dirichlet boundary condition). In addition to the hypotheses of Theorem \1.3\ assume that the 
domain C HI has boundary portion Fi which satisfies the exterior sphere condition, and that 
fGC^i^). 

(1) If 13 > 1 and also g € Cioc(ri), then the function u* in (jl.23p is a solution to problem (jl.ip 
with boundary condition ()1.3p along Ti. In particular, u* G Cioc(i^U Fi) n C^"^"(€?') and u* 
satisfies the growth assumption ()1.20p . 

(2) 7/0 < /3 < 1 and also g G C\oc{d&), then the function u* in ()1.23p is a solution to problem 
p.ip with boundary condition ()1.5p along dff. In particular, u* G Cioc(^) H C^^"(^) and 
u* satisfies the growth assumption ()1.20p . 

We next have existence of solutions to the elliptic boundary value problem when the boundary 
data g is C'^+°' up to Fi c Off. 

Theorem 1.9 (Existence of solutions to the elliptic boundary value problem with C^"*"" Dirichlet 
boundary condition). In addition to the hypotheses of Theorem \1.3l let GM be a domain such 
that the boundary portion Fi is of class (7^+", that f G Ci';^^(^U Fi) and g G Cf^+"(^U Fi). 
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(1) If P > 1, then u* , given by (jl.23p . is a solution to problem (jl.ip with boundary condition 
(jl.3p along Ti. In particular, 



u* £ c2+-(^uri) 

and satisfies the growth assumption (jl.20p . 
(2) /f < /5 < 1 and G Cioc(5^), i/ien u* , given by (jl.23p . is a solution to problem (jl.ip mt/i 
boundary condition ()1.5p along dff. In particular, 

u* G C7ioc(^)nc2+°(€?uri) 

and satisfies the growth assumption (jl.20p . 

Remark 1.10 (Existence of solutions with Daskalopoulos-Hamilton-Koch Holder regularity). When 
/ G C°(^U To) and g' = 0, we establish in |171 Theorem 1.11 & Corollary 1.13] that the solutions 
to the elhptic boundary value problem ([□]), ([L3]) lie in C'^+°'{ff U Tq) n Cioc(^) for ah (3 > 0, 
where Cg~^'^{ff UTq) is an elliptic analogue of the parabolic Daskalopoulos-Hamilton-Koch Holder 
spaces described in [HI [33]. A function u G C^"'"°(^ U Fq) has the property that u, Du,yD'^u 
are continuous up to Tq and yD^u = on Fq, where Cf{ff U Fq) is defined by analogy with 
the traditional definition of C°(^) in [25], except that Euclidean distance between points in is 
replaced by the cycloidal distance function. 

Remark 1.11 (Existence and uniqueness of solutions to elliptic boundary value problems). Existence 
and uniqueness of solutions to the elliptic boundary value problem (jl.ip and (jl.3p . provided Fi = 
dff, follow from Schauder methods when the coefficient matrix, (a*-'), of the second-order derivatives 
in A is uniformly elliptic. For example, see \25\ Theorem 6.13] for the case where is bounded 
and / and the coefficients of A are bounded and in C"(^), a G (0, 1), giving a unique solution 
u G C2+"(^) nC(^), while [25l Theorem 6.14] gives u G C2+"(#) when / and the coefficients of A 
are in C"{^). See [Ml Corollary 7.4.4], together with [Ml Corollary 7.4.9] or [Ml Theorem 7.6.4] 
or [34t Theorem 7.6.5 & Remark 7.6.6], for similar statements. 

1.2.2. Uniqueness of solutions to elliptic obstacle problems. For ^1,^2 G we set 

r-eiA6l2 



We then have the 



e~''fiZis))ds 



(1.26) 



e-^^^g(Z(ei))l|,^<,,} e-^''^^{Z{e2))l{e,<e,} 



Theorem 1.12 (Uniqueness of solutions to the elliptic obstacle problem). Let r > 0, q > 0, and f 
be as in Theorem \1.3\ and tp belong to C{ff) and satisfy (ll.20p on 0". 

(1) If 13 > I, let i/j G Cioc(^UFi) and g G Cioc(ri) obey (fT:20]) and 1^ on Fi. Let 

u G Cioc(^UFi)nC72^ 



6e a solution to the elliptic obstacle problem (jl.2p . (jl.3p stfc/i i/iat u and Au obey (jl.20p on 
^. T/ien, u = u* on UTi, where u* is given by 

n*(z) := sup j;*'^(z), (1.27) 



and Tff is defined by (jl.2ip . /or all z £ ^ U F 



1- 
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(2) IfO< p < 1, let i) e Cioc(^) and g G Qodd^) obey (fr20D and 1^ on dff. Let 

be a solution to the elliptic obstacle problem (II. 2p . (jl.Sp . such that u and Au obey (ll.20p 
on G . Then, u = u* on 6 , where u* is given by (|1.27p . 

Theorem 1.13 (Uniqueness of solutions to the elliptic obstacle problem (jl.2p . ()1.3p . when < /3 < 1). 
Letr > 0, q>0, < P <1, and f be as in Theorem\rM Let V' € Cioc(^U Ti) obey (fL20]) on G 
and let g G Cioc(ri) obey (fOOj) and ([LlD on Ti. If 

u G CxU^ U Ti) n C\G) n C]j^^{G U To) 

is a solution to the elliptic obstacle problem (11. 2p . (I1.3P such that u and Au obey (ll.20p . i/ien u = u* 
on GUTi, where u* is given by 

u*{z) := snpJ^*'''^{z), (1.28) 
and I'ff is defined by (jl.22p . /or all z £ GuTi. 

Remark 1.14 (Existence and uniqueness of strong solutions in weighted Sobolev spaces to the 
elliptic obstacle problem) . Existence and uniqueness of strong solutions in weighted Sobolev spaces 
to problem (jl.2p with Dirichlet boundary condition (|1.3p along Fi, for all /3 > 0, is proved in [9l 
Theorem 1.6], and Holder continuity of such solutions up to boundary portion Tq is proved in \15\ 
Theorem 1.13]. 

1.2.3. Existence and uniqueness of solutions to parabolic terminal/boundary value problems. We 
shall need to appeal to the following analogue of Hypothesis 11.21 

Hypothesis 1.15 (Growth condition). If f is a function then, for all {t,x,y) in its domain of 
definition, 

\v{t,x,y)\ < C7(l + e*'^i^ + e*-^2X)^ ^1_29) 

where C > 0, < A/i < ^, and M2 G [0, 1]. 

We let C{Q) denote the vector space of continuous functions on Q, while C{Q) denotes the 
Banach space of functions which are uniformly continuous and bounded on Q. We let Du denote 
the gradient and let D^u denote the Hessian matrix of a function u on Q with respect to spatial 
variables. We let C^{Q) denote the vector space of functions, u, such that u, ut, and Du are 
continuous on Q, while C^{Q) denotes the Banach space of functions, u, such that u, ut, and Du 
are uniformly continuous and bounded on Q; finally, C'^{Q) denotes the vector space of functions, 
u, such that ut, Du, and D^u are continuous Q, while C'^{Q) denotes the Banach space of functions, 
u, such that u, ut, Du, and D^u are uniformly continuous and bounded on Q. If T ^ dQ is a 
relatively open set, we let Cioc(Q U T) denote the vector space of functions, u, such that, for any 
precompact open subset F d Q U T, we have u G C{V). 

Theorem 1.16 (Uniqueness of solutions to the parabolic terminal/boundary value problem). Let 
f belong to C{Q) and obey (jl.29p . Then 

(1) > 1, assume g G Cioc(S^Q) obeys (fr29]) on S^Q. Let 



uG aoc(QugiQ)nC72(g) 
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be a solution to the parabolic terminal/boundary value problem (jl.7p . (jl.9p which obeys 
(jl.29p on Q. Then, u = u* on QU d^Q, where u* is given by 

fAT 



u*{t,z) 



t,z 



(1.30) 



-'^'-'^f{s,Z{s))ds 

Jt 

and Tff is defined by ()1.2ip . /or a// (t, z) G Q U S^Q. 
(2) 7/0 < /3 < 1, assume g e Cioc(SQ) o6eys ()1.29p on (5Q, and let 

uG aoc(QugQ)nC72(Q) 

6e a solution to the parabolic terminal/boundary value problem (II. 7p . ()1.12p which obeys 
(I1.29P on Q. Then, u = u* on QU dQ, where u* is given by ()1.30p . 

By analogy with [8l Definition 2.2], we let C^ ioc((0'^) X (^U To)) denote the subspace of 
C^'^{{0,T) X ff) nC[^^((0,T) X (i^UFo)) consisting of functions, u, such that, for any precompact 
open subset V (s (0,r) x (ffUTo), 

sup \u{t, z)\ + \ Du{t, z)\ + \y u{t, z)\ < oo. (1.31) 

Remark 1.17 (On the definition of Clf^^{{0,T) x (^U Tq)). The definition of the function space 
Cl'l^^{{0,T) X 1^ U To) can be relaxed in our paper from that implied by [HI Definition 2.2]. It 
is enough to say that a function u belongs to Cl'l^{{0,T) x ^ U Tq) n Cioc{{0,T) x ^ U Tq) if 
u G Ci'i((0, T) X ff) and (fL3T]l holds on any precompact open subset C/ G (0,r)x d ^UTq. 

We have the following alternative uniqueness result. 

Theorem 1.18 (Uniqueness of solutions to the parabolic terminal/boundary value problem when 
< /3 < 1). Let < ^ < 1 and f be as in Theorem[LM Let g G Cioc(S^Q) obey (fLMjl on ^^Q, 
and 

u G CxUQ u s'Q) n c2(Q) n c]i,((o, r) x (^ u Fq)) 

6e a solution to the parabolic terminal/boundary value problem (jl.7p . (II. 9p which obeys ()1.29p on 
Q. Then, u = u* on QU S^Q, where u* is given by 

ruffAT 1 
e-<'-'^f{s,Z{s))ds 



u*{t,z) := E; 



+ E, 



-<-^^^-^)g(u^AT,Z{i.^^AT)) 



(1.32) 



and is defined by IH^llh . for all {t, z) gQU d^Q. 

Remark 1.19 (Existence and uniqueness of strong solutions in weighted Sobolev spaces to the para- 
bolic terminal/boundary value problem). Existence and uniqueness of strong solutions in weighted 
Sobolev spaces to problem (jl.7p with Dirichlet boundary condition (II. 9p along S^Q, for all /3 > 0, 
is proved in [TO] . 



Remark 1.20 (Growth of solutions to parabolic terminal/boundary value problems). Karatzas and 
Shreve allow faster growth of solutions when the growth on the coefficients of the differential 
operator is constrained j30i Theorem 4.4.2 &: Problem 5.7.7], and polynomial growth of solutions 
is allowed for linear growth coefficients and source function / with at most polynomial growth [301 
Theorem 5.7.6]. 
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Remark 1.21 (Barrier option pricing and discontinuous terminal/boundary conditions). In appli- 
cations to finance, ^ will often be a rectangle, {xq,xi) x (0,oo), where — oo < xq < xi < oo; the 
growth exponents will be Mi = and M2 = 1 — indeed, the source function / will always be zero 
and the spatial boundary condition function g : (0, T) x Fi — )■ M will often be zero. However, the 
spatial boundary condition, g : (0, T) x Fi — )• M, and terminal condition, g : {T} x — )• M, may be 
discontinuous where they meet along {T} x dff, as in the case of the down-and-out put, with 




0<t<T,x = xo,y>0, 
t = T,xq < X < 00, y > 0, 



where g is discontinuous at {T,XQ,y) if i^T — > 0, that is, xq < log if. We shall consider 
the question of establishing stochastic representations for solutions to parabolic terminal/value 
problems (European-style option prices) or parabolic obstacle problems (American-style option 
prices) with discontinuous data elsewhere. 

For a G (0, 1), we let C'^{Q) denote the subspace of C{Q) consisting of locally a-Holder contin- 
uous functions, u, on Q, that is, for any precompact open set V ^ Q, 

r 1 \u{t\z^)-u{t^,Z^)\ 

Nc"(V) •= sup Y— 7. _ ^.r> < 00, (1.33) 



(i%^»^)ey i^\z^ - + - t^\j 



and we let C°'{Q) C C{Q) denote the Banach space of functions, u, which are uniformly Q-Holder 
continuous on Q, that is 

Mc-{Q) < 00 • 

When Q is unbounded, we let C^^{Q) denote the subspace of C"{Q) consisting of functions, u, 
such that, for any precompact open set V Q, we have 

Mc"(y) < 

We let C'^~^"{Q) denote the subspace of C^i^Q) consisting of functions, n, such that n, ut, and the 
components of Du and D^u belong to C'^{Q), and let C'^~^°'{Q) C C^^Q) denote the Banach space 
of functions, u, such that u, ut, and the components of Du and D'^u belong to C'^{Q). 

Theorem 1.22 (Existence of solutions to the parabolic terminal/boundary value problem with 
continuous Dirichlet boundary condition). In addition to the hypotheses of Theorem M.KA let C HI 
he a domain such that the boundary Fi obeys an exterior sphere condition, and f G C^^{Q). 

(1) If /3 > 1 and g € Cioc(S"'^Q)> then u* in (jl.30p is a solution to problem (jl.7p with bound- 
ary condition (jl.9p . In particular, u* G Cioc(Q U d^Q) H C'^~^°'{Q) and obeys the growth 
assumption (jl.29p . 

(2) /f < /? < 1 and g € Cioc(SQ)) then u* in ()1.30p is a solution to problem (jl.7p with 
boundary condition (jl.l2p . In particular, u* £ Cioc(Q U (5Q) Ci C^^"((5) and satisfies the 
growth assumption ()1.29p . 

For T ^ (5Q a relatively open subset, we let Cf+"(Q U T) denote the subspace of C^+'^{Q) such 
that, for any precompact open set U (£ Q UT, we have u G C'^'^'^ (tj) . 

Theorem 1.23 (Existence of solutions to the parabolic terminal/boundary value problem with 
(j2+a j])ij-ic]2iet boundary condition). In addition to the hypotheses of Theorem \1.16[ let iff cM. be 
a domain such that the boundary portion Fi is of class C^^". 
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(1) ///3 > 1 andge Cf+°(Q U S^Q) obeys 

-gt + Ag = f on{T}>cVi. (1.34) 
Then u* in (ll.30p is a solution to problem (jl.7p with boundary condition (jl.9p . /n particular, 

u* eClt^iQuQ'Q) 
and obeys the growth estimate (ll.29p . 

12+c 
loc 

- 5t + ^5 = / on {T} X a^. (1.35) 
T/ien u* in (jl.30p is a solution to problem ()1.7p wii/i boundary condition ()1.12p . /n partic- 



{2) If0</3<1 andg£ Cl+''{Q U S^Q) n Cioc(Q) oteys 



ular, 



u* GCl+''{QU(5'Q)nCioc{Q)- 



""loc 

and o6eys i/ie growth estimate (|1.29p . 

Remark 1.24 (Zero and first-order compatibility conditions for parabolic equations). The conditions 
(|1.34p and (|1.35p are the analogues of the first-order compatibility condition |34| Equation (10.4.3)]. 
The analogue of the zero-order compatibility condition in j34( Equation (10.4.2)] automatically 
holds at {T} x Fi or {T} x dff, since we always choose h = g{T, ■) on Fi or dff, respectively, in 
this article. 

Remark 1.25 (Existence of solutions with Daskalopoulos-Hamilton-Koch Holder regularity). As in 
the elliptic case, the solutions to the parabolic terminal/boundary value problem (11. 7p . (11. 9p should 
lie in Cf+"(Q) n CiodQ) for all /3 > 0, where C^+"(Q) is the parabohc Daskalopoulos-Hamilton- 
Koch Holder space described in [TT1[33]. A function u G Cl^°'{Q) has the property that n, Du, yD^u 
are continuous up to Fq and yD^u = on (0, T) x Fq, where C^{Q) is defined by analogy with 
the traditional definition [34] of C°'{Q), except that Euclidean distance between points in Q is 
replaced by the cycloidal distance function. When Q = H x (0, T), we establish this existence 
result in [18^ Theorem 1.1]. 

Remark 1.26 (Existence and uniqueness of solutions to parabolic terminal/boundary value prob- 
lems). Existence and uniqueness of solutions to the parabolic terminal/boundary value problem 
(II. 7p and (II. 9p . again provided Fi = dff, follow from Schauder methods when the coefficient matrix, 
(a*-'), of A is strictly elliptic on ff. For example, see [351 Theorems 5.9 & 5.10] for the case where / 
and the coefficients of A are bounded and in C"((5), giving a unique solution u € C^~'"°(Q) r\C{Q). 

1.2.4. Uniqueness of solutions to parabolic obstacle problems. For 01,02 G ^t,T, < t < T, we set 



-^^'-'^f{s,Z{s))ds 



('^~'^{d2,Z{02))l{e,<e,} 



(1.36) 



e-'-(^-*)<7(^i,^(^i))l{e,<e.} 



We have the following uniqueness result of solutions to the parabolic obstacle problem with different 
possible boundary conditions, depending on the value of the parameter /? > 0. 

Theorem 1.27 (Uniqueness of solutions to the parabolic obstacle problem). Let f be as in Theorem 
\1.16\ and ip belong to C{Q) and satisfy (jl.29p . 

(1) If P > I, assume tp G Cioc((5 U S^Q) and g G Cioc(S^Q) obeys ^L29\i on d^Q and (fTTTD . 
Let 

«GCioc(QugiQ)nc2(g) 
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be a solution to the parabolic obstacle problem (jl.Sp . (jl.9p such that u and Au obey ()1.29p 
on Q. Then, u = u* on QuS^Q, where u* is given by 

u*{t,z):= sup J^'^^^'\t,z), (1.37) 

and Tff is defined by (jl.2ip . /or all (t, z) £ Q U d^Q. 
(2) IfO< p < 1, assume ^ G CiodQ) and g £ CiodSQ) obeys I^L29\\ on SQ and (fTTiD . iet 

uG C7i„,(QugQ)nC72(Q) 

6e a solution to the parabolic obstacle problem (jl.Sp . (jl.l2p suc/i i/iai u and Au obey (|1.29p 
on Q. Then, u = u* on QU dQ, where u* is given by (|1.37p . 

Theorem 1.28 (Uniqueness of solutions to the parabolic obstacle problem (jl.Sp . (jl.9p . when 
< /? < 1). Lei < /? < 1 and / be as in Theorem \1.16[ Assume ip G Cioc(Q U ^^Q)j and 
g G Cioc(3^(5) obey (fL29D on S^Q and (fLTT]l . iei 

G Cioc(Q u giQ) n c2(Q) n cI'Uq u (o, t) x (^ u Fq)) 

be a solution to the parabolic obstacle problem (jl.Sp . ()1.9p such that u and Au obey ()1.29p . Then, 
u = u* on QU d^Q, where u* is given by 

u*{t,z):= sup J^<'^^'^(i,z), (1.38) 
and Vff is defined by (ll.22p . /or a// (t, z) G Q U S^Q. 

Remark 1.29 (Existence and uniqueness of strong solutions in weighted Sobolev spaces to the 
parabolic obstacle problem). Existence and uniqueness of strong solutions in weighted Sobolev 
spaces to problem (II. 8p with Dirichlet boundary condition (II. 9p along S^Q, for all /3 > 0, is proved 
in ^U\. 

1.3. Survey of previous results on stochastic representations of solutions to boundary 
value or obstacle problems. Stochastic representations of solutions to elliptic and parabolic 
boundary value and obstacle problems discussed by Bensoussan and Lions [5] and Friedman [23] 
are established under the hypotheses that the matrix of coefficients, (a*-'), of the second-order 
spatial derivatives in an elliptic linear, second-order differential operator, A, is strictly elliptic and 
that all coefficients of A are bounded. Relaxations of these hypotheses, as in [241 Chapter 13 &: 15], 
and more recently [47], fail to include the Heston generator mainly because the matrix {a^-') does 
not satisfy 

Hypothesis 1.30 (Extension property for positive definite, matrix- valued functions). Given a 
subdomain V ^ (0,oo) x M'^, for d> 1, we say that a matrix-valued function, 

a:V ^ M'^^'^, 

which is on V and a{t, z) is positive definite for each (t, z) £ V has the extension property if 
there is a matrix-valued function, 

d : [0,oo) xR'^ ^ R'^'"^, 

which coincides with a on V but is on [0, oo) X M and d{t, z) is positive definite for each 
(t,z) G [0,oo) X M'^. 
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Naturally, Hypothesis 11.301 is also applicable when the matrix a is constant with respect to time, 
that is, in elliptic problems. Note that in the case of the Heston process, d = 2, V = {0, oo) x H, 
and 



and so the matrix a does not satisfy Hypothesis 11.301 We now give more detailed comparisons for 
each of the four main problems which we consider in this article. 

1.3.1. Elliptic boundary value problems. Uniqueness of stochastic representations of solutions to 
non- degenerate elliptic partial differential equations is established in [24^ Theorem 6.5.1], |30l 
Proposition 5.7.2], [Ml Theorem 9.1.1 & Corollary 9.1.2, Theorem 9.3.2], and [H Theorem 2.7.1 & 
Remarks 2.7.1, 2.7.2] (for a bounded domain ff), and (Sj Theorem 2.7.2 Sz Remarks 2.7.3-5] (when 
the domain is the whole space, = M"). 

Existence (and uniqueness) of stochastic representations of solutions to non- degenerate elliptic 
partial differential equations is established in [241 Theorem 6.5.1], [381 Theorem 9.2.14], and [381 
Theorem 9.3.3 &: Remark, p. 196]. 

Existence and uniqueness of solutions to a certain class of degenerate elliptic partial differential 
equations are described by Friedman in |24l Theorems 13.1.1 &: 13.3.1], but those results do not 
apply to the Heston operator because a square root, (o"*-'), of the matrix (a*-') cannot be extended as 
a uniformly Lipschitz continuous function on M^, that is, \2A\ Condition (A), p. 308] is not satisfied. 
Stroock and Varadhan j33j §5-8] also discuss existence and uniqueness of solutions to degenerate 
elliptic partial differential equations, but their assumption that the matrix (a*-') satisfies Hypothesis 
11.301 does not hold for the Heston operator (see \4:3\ Theorem 2.1]). 

More recently, Zhou [47] employs the method of quasiderivatives to establish the stochastic 
representation of solutions to a certain class of degenerate elliptic partial differential equations, 
and obtains estimates for the derivatives of their solutions. However, his results do not apply 
to the Heston operator because [47^ Assumptions 3.1 & Condition (3.2)] are not satisfied in this 
case. Moreover, the Dirichlet condition is imposed on the whole boundary of the domain (see 
[171 Equation (1.1)]), while we take into consideration the portion of the boundary, Tq, where the 
differential operator A becomes degenerate. 

1.3.2. Elliptic obstacle problems. We may compare Theorems 11.121 and 11.131 with the uniqueness 
assertions (in increasing degrees of generality) for non- degenerate elliptic operators in [5l Theorems 
3.3.1, 3.3.2, 3.3.4, 3.3.5, 3.3.8, 3.3.19, 3.3.20, & 3.3.23]. See also [Ml Theorem 10.4.1] and [Ml The- 
orems 16.4.1, 16.4.2, 16.7.1, &: 16.8.1] for uniqueness assertions non- degenerate elliptic operators, 
though with more limited applicability. 

1.3.3. Parabolic boundary value problems. Uniqueness of solutions to non-degenerate parabolic par- 
tial differential equations and their stochastic representations are described in [241 Theorems 6.5.2, 
6.5.3], [Ml Theorem 5.7.6] and [Sj Theorems 2.7.3 & 2.7.4]. 

Friedman obtains fundamental solutions and stochastic representations of solutions to certain 
degenerate parabolic partial differential equations in |23] , while he obtains uniqueness and stochastic 
representations of solutions to the Cauchy problem in |22] : those results are summarized in [241 
Chapter 15]. Nevertheless, the results in [241 Chapter 15] and [23] do not apply to the Heston 
operator because Hypothesis 11.301 does not hold, that is [24l Condition (A), p. 389] is not satisfied. 
Therefore, the method of construction in [231 Theorem 1.2] of a candidate for a fundamental solution 
does not apply to the Heston operator. A stochastic representation for a solution to the Cauchy 
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problem for a degenerate operator is obtained in |24l §15.10], but the hypotheses of [24:\ Theorem 
15.10.1] are again too restrictive and exclude the Heston operator. 

Ekstrom and Tysk |14j consider the problem of pricing European-style options on an underly- 
ing process which is the solution to a degenerate, one-dimensional stochastic differential equation 
which satisfies Hypothesis 2.1], and so includes the Feller square root (or Cox-Ingersoll-Ross) 
process, (12. ip . The option price is the classical solution in the sense of |14^ Definition 2.2] to the 
corresponding parabolic partial differential equation [141 Theorem 2.3]. Under their assumption 
that the payoff function g{T,-) is in C^([0, oo)), they show that their classical solution has the 
regularity property, 

ue c([o,r] X [o,oo))nc^([o,r) x [0,00)) nC72([o,r) x (0,00)), 

and obeys the second-order boundary condition, 

lim yuyy{t,y)=0, VioG(0,r) (by [Ml Proposition 4.1]). 

As a consequence, in the framework of our article, their solution obeys 

^ e Cl'Lm to) X [0, oo)), Vto G (0, T), 

where the vector space of functions C'^'ioc((^' ^o) x [0, oo)) is defined by analogy with ()1.3ip . 

In |13] . Ekstrom and Tysk extend their results in [14] to the case of two-dimensional stochastic 
volatility models for option prices, where the variance process satisfies the assumptions of [141 
Hypothesis 2.1]. 

Bayraktar, Kardaras, and Xing [3] address the problem of uniqueness of classical solutions, in 
the sense of [H Definitions 2.4 & 2.5], to a class of two-dimensional, degenerate parabolic partial 
differential equations. Their differential operator has a degeneracy which is similar to that of 
the Heston generator, —A, and to the differential operator considered in [T3], but the matrix of 
coefficients, (a*-'), of their operator may have more than quadratic growth with respect to the spatial 
variables (see [H Standing Assumption 2.1]). Therefore, weak maximum principles for parabolic 
partial differential operators on unbounded domains such as [Ml Exercise 8.1.22] do not guarantee 
uniqueness of solutions in such situations. The main result of their article - [H Theorem 2.9] - 
establishes by probabilistic methods that uniqueness of classical solutions, obeying a natural growth 
condition, holds if and only if the asset price process is a martingale. 

In our article, we consider the two-dimensional Heston stochastic process, (|1.17p . where the 
component Y of the process satisfies [141 Hypothesis 2.1] and [H Standing Assumption 2.1]. We 
only require the payoff function, g{T,-), to be continuous with respect to the spatial variables 
and have exponential growth, as in ()1.29p . Notice that the conditions on the payoff function are 
more restrictive in [141 Hypothesis 2.1] and [H Standing Assumption 2.3] than in our article. We 
consider the parabolic equation associated to the Heston generator, —A, on bounded or unbounded 
subdomains, ff, of the upper half plane, H, with Dirichlet boundary condition along the portion, 
Fi, of the boundary dff contained in H. Along the portion, Tq, of the boundary contained in 
dM, we impose a suitable Dirichlet boundary condition, depending on the value of the parameter 
/3 in (ll.lSp . which governs the behavior of the Feller square-root process when it approaches the 
boundary point y = 0. In each case, we establish uniqueness of solutions by proving that suitably 
regular solutions must have the stochastic representations in Theorems 11.161 and 11.181 and we prove 
existence and regularity of solutions, in a special case, in Theorems 11.221 and 11.231 complementing 
the results of [14] . In addition, we consider the parabolic obstacle problem and establish uniqueness 
and the stochastic representations of suitably regular solutions in Theorems 11.271 and 11.281 
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1.3.4. Parabolic obstacle problems. We may compare Theorems 11.271 and 11.281 with the uniqueness 
assertions and stochastic representations of solutions (in increasing degrees of generality) for non- 
degenerate operators in [5l Theorems 3.4.1, 3.4.2, 3.4.3, 3.4.5, 3.4.6, 3.4.7, 3.4.8]. 

1.4. Further work. The authors are developing an extension of the main results of this article to 
a broader class of degenerate Markov processes in higher dimensions and more general boundary 
conditions (including Neumann and oblique boundary conditions). 

1.5. Outline of the article. For the convenience of the reader, we provide a brief outline of 
the article. We begin in f|2] by reviewing or proving some of the key properties of the Feller 
square root and Heston processes which we shall need in this article. In ^ we prove existence 
and uniqueness (in various settings) of solutions to the elliptic boundary value problem for the 
Heston operator, while in ^ we prove uniqueness (again in various settings) of solutions to the 
corresponding obstacle problem. We proceed in ^ to prove existence and uniqueness of solutions 
to the parabolic terminal/boundary value problem for the Heston operator and in ^ we prove 
uniqueness of solutions to the corresponding parabolic obstacle problem. Appendices |A] and [B] 
contain additional technical results which we shall need throughout our article. 

1.6. Notation and conventions. When we label a condition an Assumption, then it is considered 
to be universal and in effect throughout this article and so not referenced explicitly in theorem and 
similar statements; when we label a condition a Hypothesis, then it is only considered to be in 
effect when explicitly referenced. We let N := {1, 2, 3, . . .} denote the set of positive integers. For 
X, y G M, we denote x Ay := min{x, y}, xM y := max{x, y} and x~^ := xM 0. 

1.7. Acknowledgments. We are very grateful to everyone who has provided us with comments 
on previous versions of this article or related conference or seminar presentations. Camelia Pop 
thanks Daniel Ocone for many helpful discussions on probability theory. Finally, we thank the 
anonymous referee for many helpful suggestions and kind comments. 

2. Properties of the Heston stochastic volatility process 

In this section, we review or develop some important properties of the Feller square root process 
and the Heston stochastic volatility process. 

By [161 Theorem 1.9], it follows that for any initial point {t,y) G [0, oo) x [0,oo), the Feller 
stochastic differential equation, 

dY{s) = k{^- Y{s)) ds + a^J\Y{s)\dW{s), s > t, 
Yit) =y, 

admits a unique weak solution (Y*'''^ (s),W{s))s>t, called the Feller square root process, where 
{W{s))s>t is a one-dimensional Brownian motion on a filtered probability space (Q, ^,P*'^,F) such 
that the filtration F = {^(s)}s>o satisfies the usual conditions |30l Definition 1.2.25]. Theorem 
1.9 in [16] also implies that the Heston stochastic differential equation ()1.17p admits a unique weak 
solution, (Z*'^(s), Ty(s)) for any initial point {t,z) G [0, cxd) x H, where {W{s))s>t is now an 

M^-valued Brownian motion on a filtered probability space (J7, Q*'^, F) such that the filtration 
F = {^{s)}s>o satisfies the usual conditions. When the initial condition {t, y) or {t, z) is clear from 
the context, we omit the superscripts in the definition of the probability measures P*'^ and Q*'^, 
respectively. 

Moreover, the weak solutions to the Feller and Heston stochastic differential equations are strong. 
To prove this, we begin by reviewing a result of Yamada |46j . 
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Definition 2.1 (Coefficients for a non-Lipscliitz stocliastic differential equation). [461 p. 115] In 
this article we shall consider one-dimensional stochastic differential equations whose diffusion and 
drift coefficients, a, b, obey the following properties: 

(1) The functions a, b : [0, oo) x R — t- M are continuous. 

(2) (Yamada condition) There is an increasing function g : [0, oo) — >■ [0, oo) such that q{0) = 0, 
for some e > one has Q~'^{y) dy = oo, and 

\a{t,yi)-a{t,y2)\<Q{\yi-y2\), yi,y2GM,t>0. (2.2) 

(3) There is a constant Ci > such that 

\h{t,y2)-b{t,yi)\<Ci\y2-yil yi,y2GM,t>0. (2.3) 

(4) There is a constant C2 > such that 

\a{t,y)\ + \b{t,y)\<C2{l + \y\), t > 0, y G M. (2.4) 

Clearly, the coefficients of the Feller stochastic differential equation obey the hypotheses in 
Definition 12.31 where a{t,y) = Oy/y and b{t,y) = k{9 — y). Indeed, one can choose Ci = k, 
C2 = max{K, k9, ct}, and Q{y) = cjy/y., as the mean value theorem yields 

Vm- Vvi = c{yi,y2){y2 - yi), 

where 

N I I 1 

c{yi,y2) = 7; I , =^ < -r 

2^0 y^yi+t{y2-yi) Vy2 - yi 

for < yi < y2- See j46[ Remark 1] for other examples of suitable functions g. 

Remark 2.2. When g{u) = u', 7 E [^!l] [461 Remark 1], then Definition 12.11 implies that a(t, •) is 
Holder continuous with exponent 7, uniformly with respect to t S [0,oo). 

Definition 2.3 (Solution to a non-Lipschitz stochastic differential equation). |46[ p. 115], [401 
Definitions IX. 1.2 & IX. 1.5] Let (r2,^,P, F) be a filtered probability space satisfying the usual 
conditions. We call a pair {Y {s) ,W {s)) s>q a weak solution to the non-Lipschitz one-dimensional 
stochastic differential equation, 

dY{s) = b{s,Y{s))ds + a{s,Y{s))dW{s), s>0,y(0)=y, (2.5) 

where y € M, if the following hold: 

(1) The processes Y{s) and W{s) are defined on (r2,^,P, F); 

(2) The process Y{s) is continuous with respect to s S [0, 00) and is F-adapted; 

(3) The process W{s) is a standard F-Brownian motion. 

We call {Y{s), W{s))s>q a strong solution to ()2.5p if Y is F'^-adapted, where F^ is the P-completion 
of the filtration of ^ generated by (VF(s))s>o. (Compare j28[ Definition IV. 1.2], [301 Definition 
5.2.1], and [Ml §5.3].) 

Theorem 2.4. [l6l p. 117] There exists a weak solution {Y,W) to ()2.5p . 

Remark 2.5. Yamada's main theorem [461 p. 117] asserts considerably more than Theorem 12. 4[ 
In particular, his article shows that ()2.5p may be solved using the method of finite differences. 
Simpler results may suffice to merely guarantee the existence of a weak solution, as we need here; 
see Skorokhod [42] . 

Proposition 2.6. There exists a unique strong solution to (|2.5p . 
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Proof. Theorem 12.41 ensures that (j2.5|) admits a weak solution. Conditions (j2.2p and (j2.3p ensure 
that pathwise uniqueness holds for (weak) solutions to (j2.5p by Revuz and Yor |40l Theorem IX. 3. 5 
(ii)], while Karatzas and Shreve \30\ Corollary 5.3.23] imply that (j2.5p admits a strong solution; 
see [30l p. 310]. Conditions (12. 2 p and ()2.3p guarantee the uniqueness of strong solutions to p.Sp by 
Karatzas and Shreve [3^ Proposition 5.2.13]; compare Yamada and Watanabe (Pathwise 
uniqueness is also asserted for (12. 5 p by \28\ Theorem IV. 3. 2] when (12. 5 p is time-homogeneous, noting 
that the coefficients a, b are not required to be bounded by Ikeda and Watanabe [28\ p. 168]). We 
conclude that a strong solution to (12. 5p exists and is unique. □ 



Corollary 2.7. Given any initial point {t, y) G [0, oo) x [0, oo), there exists a unique strong solution, 
(Y*'y{s),W{s))s>t, to the Feller stochastic differential equation. 

Proof. Immediate from Proposition 12.61 □ 

Corollary 2.8. Given {t,z) G [0, oo) xH, there exists a unique strong solution, (Z^''^{s),W{s))s>t, 
to the Heston stochastic differential equation, where {W{s))s>() is a standard two-dimensional F- 
Brownian motion on (r2,^,P,F). 

Proof. By Proposition 12.61 the Cox-Ingersoll-Ross stochastic differential equation has a unique 
strong solution, (y*'^(s), VF2(s))s>t) where (W2{s))s>t is a standard one-dimensional F-Brownian 
motion on (f],^,P,F) and {Y^'yls))s>t is F^^-adapted. But given {Y^^y{s))s>t and a standard 
two-dimensional F-Brownian motion, {W{s))s>t = iWi{s),W2{s))s>t on [Q, ^ ,¥,¥), the process 
(X^'^'y{s))s>t, and thus {Z^'%s))s>t = (X*'^'2^(s), y*'2'(s)),>t, is uniquely determined by 



X*'^'?/(s) =x + (^r-q- ^y*'^(n)^ 



du 



+ 



This completes the proof. □ 

Lemma 2.9 (Properties of the Feller square-root process). The unique strong solution of the Feller 
stochastic differential equation started at any {t, y) G [0, oo) x [0, oo) satisfies 

Y{s) > P*'^-a.s., Vs > t, (2.6) 

and also 



i: 



l{y(„)=o} du = 0, Vs > t, (2.7) 
L(s, x) =0, Vx < 0, Vs > t, (2.8) 
where L(-, ■) is the local time of the Feller square-root process. 

Proof. Without loss of generality, we may assume that t = 0. In [H Lemma 2.4], it is proved that 
L(s, 0) = 0, for all s > 0, but it is not clear to us why it also follows that 

L(s, 0-) := lim L(s, x) = 0, Vs > 0, 

a property we shall need in our proof of (12. 6p . To complete the argument, we consider the following 
stochastic differential equation, 

dY{s) = b{Y{s)) ds + a{Y{s)) dW{s), s > 0, 

Y{0)=y, 



18 



P. FEEHAN AND C. POP 



where we let 

b{y) := - y) and a{y) := l{y>o}a^, Vy G M. (2.9) 

This stochastic differential equation admits a unique strong solution by Proposition I2.6i We will 
show that Y{s) > a.s., for all s > 0, so that uniqueness of solutions to the Feller stochastic 
differential equation (12. ip implies that Y = Y a.s. and Y will satisfy the same properties as Y . 
Thus, it is enough to prove (j2.7p and (j2.8p for Y . Property (j2.6p is a consequence of the preceding 
two properties of Y . 

Let L be the local time process for the continuous semimartingale Y (see [301 Theorem 3.7.1]). 
Prom [301 Theorem 3.7.1 (iii)], we know that, for any Borel measurable function A; : M — > [0,oo), 
we have 

k{Y{u))a^Y+{u)du = 2 k{x)L{s,x)dx, Vs > 0. (2.10) 



JR 
Assume, to obtain a contradiction, that L(s, 0) > 0. Prom the right-continuity in the spatial variable 
of L{s,-) \30\ Theorem 3.7.1 (iv)], there are positive constants c and xq such that L{s,x) > c, for 
all X € [0,xo]. Por e > 0, we define k{x) = x~^, for x € [e, xq], and otherwise. With this choice 
of k, the left-hand-side in identity ()2.10p is bounded in absolute value by cj^s, for any e > 0, while 
the right-hand-side of ()2.10p is greater or equal than 2c log (xq/e), which diverges as e tends to 0. 
Therefore, our assumption that L{s, 0) > is false, and so L(s, 0) = 0. Moreover, we notice that for 
any bounded, Borel-measurable function k with support in (— oo,0) the left-hand-side in identity 
()2.10p is identically zero. Thus, we conclude that L{s,x) = 0, for all x < 0, and also L(s,0— ) = 0. 

We use this result to show that ¥(Y[s) < 0,Vs > 0) = 0. Prom [30 [ p. 223, third formula] and 
the fact that k, t9 > 0, we see that 

s 



= L{s,0) - L{s,0-) = l^y^^^^^^du 



which implies that F{Y{s) = 0, Vs > 0) = 0. It remains to show that F(Y{s) € (-oo, 0)) = 0, for all 
s > 0, which is equivalent to proving that for any e > and s > 0, we have P(y(s) G (— oo, — e)) = 0. 
Let : M — 7- [0, 1] be a smooth cut-off function such that 99|(_oo^_£) = 1 and 9j|(o,oo) = 0- We can 
choose (f such that if' < 0. Then, it follows by Ito's formula that 

viY{s)) = v?(y(0)) + r (^i^ - Y{u))^'iYiu)) + la\Y{uW{Y{u))] du 







2 

+ / a{Y{u))ip'{Y{u))dW{u) 



ip{Y{0))+ -Y{u))ip'{Y{u))du {as a{y) = when (p' ^ 0) . 





We notice that the right-hand-side is non-negative, while the left-hand-side is non-positive, as 
99' < on M, and ip' = on (0, 00). Therefore, we must have ip{Y{s)) = a.s. which implies that 
¥{Y{s) G (—00, —e)) = 0. This concludes the proof of the lemma. □ 

For a,y,t > 0, we let 

Tj^'y := inf {s>t: y*'J'(s) = a) (2.11) 

denote the first time the process Y started at y at time t hits a. When the initial condition, {t,y), 
is clear from the context, we omit the superscripts in the preceding definition (|2.1ip . Also, when 
t = 0, we omit the superscript t. 
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Lemma 2.10 (Boundary classification at ?/ = of the Feller square root process). Let he 

the unique strong solution to the Feller stochastic differential equation (|2.ip with initial condition 
yy(0) = y. Then 

(1) For 13 > 1, y = is an entrance boundary point in the sense of 132\ §15. 6(c)]. 

(2) For 0</3<l,y = 0isa regular, instantaneously reflecting boundary point in the sense of 
[321 §15.6(a)], and 

limr(5' = a.s., (2.12) 

yiO 

where Tq is given by (|2.1ip . 

Proof A direct calculation give us that the scale function, s, and the speed measure, m, of the 
Feller square root process are given by 

2 

where /3 = iK'^ja^ and /U = iKja^ . We consider the following quantities, for < a < 5 < oo and 
a; > 0, 



s(y) = y-f^e^y and m(y) = — /-^e-'^S', Vy > 



S'[a,6]:= / 5{y)dy, S {a, b] := lira S[c,b], 

J a c4.a 



M[a,6]:= / m{y)dy, M(a, 6] := lim M[c, 6], 

N{0) := r S[y,x]m{y)dy. 
Jo 

Then, for /3 > 1, we have 5(0,2;] = oo and A^(0) < oo, which implies that y = is an entrance 
boundary point ((321 p. 235]), while for < /3 < 1, we have 5(0, x] < oo and M{0,x] < oo, and so 
y = is a regular boundary point ( j321 p. 232]). 

Next, we consider the case < /3 < 1. To establish (|2.12p . we consider the following quantities 

:=r"m<T.) = fM. 

f-b ry 

Va,b{y) ■= [Ta A Tfo] = 2ua,b{y) / S[z, b]m{z)dz + 2(1- Ua,b{y)) I S{a, z]m{z)dz, 

Jy J a 

as in [321 Equations (15.6.1) & (15.6.5)] and [32l Equations (15.6.2) & (15.6.6)], respectively. Notice 
that Ta — T- Tq , when y J, 0, by the continuity of the paths of Y. Then, for fixed 6 > 0, we obtain 

limP^(r6 < To) = limlimPJ'(Tb < Ta) = 0, 

ylO ylO aiO 

hm El [To A Tb] = lim hm [Ta A Tb] = 0, 

yiO yiO alO 

from where (12. 12^ follows. □ 

Next, we have the following 

Lemma 2.11 (Properties of the Heston process). Let (Z(s))^>q be the unique strong solution to 
the Heston stochastic differential equation (I1.17|) . 

(1) Assume q >0 and r G M. Then, for any constant c G [0, 1], 

-rcs^cX{s)\ is a positive supermarting ale. (2.13) 
J s>Q 



e 
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(2) For any positive constant c < fi, 

/ ^-cKi9s^cY(s)\ is a positive supermarting ale. (2-14) 
V / s>0 

Proof. To establish (j2.13p . we use Ito's formula to give 

d (e-^'^'e'^^'A = -e-'^'e"^^'^ (cq + ^c(l - c)Y{s)) ds 
^ ^ V 2 / (2.15) 

+ ce-'^'e''^^''^,/Y{s)dWi{s). 

Notice that the drift coefficient is non-positive, since Y{s) > a.s. for all s > by Lemma [2.91 and 
g> 0, and cG [0,1]. 

Similarly, to establish (|2.14p for the Feller square root process, we have 

d L-^^^'e^^^'A = e-^'^'^'e'^^^'^c [ca'^ /2 - k) Y{s)ds 

^ ' , , , (2.16) 

+ me-''''^'e''^^'^^[Y{^(pdWx{s) + ^/Y^dW^is 



When c < p, we see that the drift coefficient in the preceding stochastic differential equation is 
non-negative. 

The supermartingale properties (j2.13p and (j2.14p follow if we show in addition that the processes 
are integrable random variables for each time s > 0. For simplicity, we let Q{s) denote either one of 
the processes we consider, and we let 6'„ be the first exit time of the Heston process (X(s), y(s))^>Q 
from the rectangle {—n,n) x (— n,n), where n € N. We set Qn.(s) := Q{s A On), for all s > 0. We 
then have 

dQn{s) = l{s<e^}dQn{s), Vs > 0, Vn G N. 

Using equations (|2.15p and (I2.16|) . it is clear that {Qnis))s>o are supermartingales, because the 
coefficients of the stochastic differential equations are bounded and the drift terms are non-positive. 
Therefore, we know that 

[Qn{t)l^{s)] < Qn{s), yt >s, Vs > 0, Vn G N. (2.17) 

Clearly, we also have Qn{t) — Q{t) a.s., as n — t- oo, for all t > s and s > 0. Taking the limit as 
n — >• oo in (|2.17p and using the positivity of the processes, Fatou's lemma yields 

E^'' [Q{tW{s)] < hminf E^^ [Q„(t)|^(s)] 

<liminfQ„(s) (by dim) ) 
= Q(s), Vt > s, Vs > 0, 
and so ([233]) and (l23iD follow. □ 

The next lemma is used to show that the functions u* given by ()1.23p and (jl.25p are well-defined 
and satisfy the growth assumption (jl.20p . 

Lemma 2.12. Suppose r > 0, and f, g, ip are Borel measurable functions on G and satisfy 
assumption (jl.20p . Then there is a positive constant C, depending on r, k, "d, Mi, M2 and C in 
()1.20p . such that for any 61,62 E the function jf^'^^ in (jl.26p satisfies the growth assumption, 

I Jf"''(^,y)l <c{i + e/'^y + e*^2") , yix,y) G G, 

where < Adi < min {r/ (k??) , ^u} and M2 G [0, 1) are as in (|1.20p . 

Remark 2.13. The obstacle function tp in ()1.26p is only relevant for solutions to problem (jl.2p . 
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Proof. The conclusion is a consequence of the properties of the Heston process given in Lemma 
12.111 We first estimate the integral term in (jl.26p . For z £ ff, then 



0iA6»2 



e-^ViZ{s))\ds 







< CE, 

<cli + 



(by ^M) 







^-M^K'ds^MiY{s) 



ds 







ds 



Using the condition Mi < min {r/(Ki?), ^} and (|2.14p . together with M2 < 1 and (j2.13p . we see 
that 



e-''ViZis))\ds 



<C{1 + e^'^y + e^2") 



(2.18) 



for a positive constant C depending on r, M\K'd, M2 and the constant C in the growth assumption 
(fOOl) on /, g and ^/;. 

Next, we show that the first non- integral term in (jl.26p can be written as 



-r6»i 



(2.19) 



for any 9i £ ^ which is not necessarily finite. This is reasonable because by rewriting 



EJ 



e~^'^g{Z{9,))l^e,<9,} 



+ E, 



e-^'^g{Z{ei))l{T<9,<e,} 



we shall see that the second term converges to zero, as T — > 00. Using the growth assumption on 
g in p.20p . we have 



E^ 



e~-'^\g{Z{9r))\l{T<e^<e,} 



^{T<ei} 



and so by Lemma I2.1H we obtain 



E^ 



e~'^'^g{Z{d^))l{T^e,<e,} 



< C fe"""^ + e-ir-MiKi9)T^Ahy _^ ^-r{l-M2)T ^M2x\ ^ 



Since Mi < r/(«;i9) and M2 < 1, we see that the right hand side converges to 0, as T — )■ 00. This 
justifies the identity (j2.19p . 

Now, we use Fatou's lemma to obtain the bound (jl.20p on the first non-integral term in (jl.26p . 
For z £ ff, 



E^ 



e-^'-\g{Z{ei))\l{e,<e2} 



< liminfEl 



< liminf C ( 1 + E^ 

n—>-co 



g-r(6liAn)gAfiy(6»iAn) 



+ E^ 



g-r(eiAn)gM2X(eiAn) 



(by dTlQD). 
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Because Mi < /x, we may apply the supermartingale property (j2.14p with c := M2. We use also 
that Ml < r I to obtain Miki) < r, and so it follows by the Optional Sampling Theorem [301 
Theorem 1.3.22] that 



El 



g-r(0iAn)gAfiy(9iAn) 



^-Mi Ki?(6»i An) gA/i y {6»i An) 



<E| 
< e^'^^y, Vn G N. 

Using the fact that M2 < 1, we see by the supermartingale property (j2.13p applies with c := Mi. 
By the Optional Sampling Theorem [301 Theorem 1.3.22] we have 



Therefore, we obtain 



g-r(0iAn)gM2X(0iAn) 



< El 



g-rA/2{eiAn)gAf2X(6»iAn) 



< e 



Af2a; 



Vn G N. 



< C (1 + e^'^y + e^'^^) 



We obtain the same bound on the second non-integral term in ()1.26p because the obstacle function 
satisfies the same growth condition (jl.20p as the boundary data g. □ 

To prove Theorems 11.161 and 11.181 we make use of the following auxiliary result 

Lemma 2.14. Let z G H and T £ (0,To], where Tq is a positive constant. Let {Z^{s))s>o the 
unique strong solution to the Heston stochastic differential equation (|1.17p with initial condition 
Z^(0) = z. Then there is a positive constant c, depending on y, k, 1), a and Tq, such that for any 
constant p satisfying 



0<p< 



we have 



sup E| 



2aT 

pX^0) 



< 00, 



(2.20) 
(2.21) 



where .%^t denotes the set of {ft, ,¥)- stopping times with values in [0,T]. 

Proof. We use the method of time-change. Denote 

M,{t):= [ ^/Y{^dWi{s), i = l,2, 



and observe that there is a two-dimensional Brownian motion (i?i,i?2) [301 Theorem 3.4.13] such 
that ^ 

Mi{t) = Bi(^j^Y{s)ds^ , i = l,2. 
Thus, we may rewrite the solution of the Heston stochastic differential equation (|1.17p in the form 

X{t) =x + {r-q)s-^ Y{s)ds + Bi (^J^ Y{s)ds^ , (2.22) 



Y{t) = y + K^s-K Y{s)ds + aBs Y 



{s)ds , 



(2.23) 



where -63 := pBi + y^l — p^B2 is a one-dimensional Brownian motion. 
For any continuous stochastic process {P(t))t>o, we let 

Mp{t) := max P{s), Vt > 0. 

0<s<t 
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We first prove the following estimate. 

Claim 2.15. There are positive constants uq and c, depending on y, k, 'd, a and Tq, such that 

QMn<My(T)<n + l)<4=e"'"/^'"'''^l{n,>no} + l{n<no}, VnGN. (2.24) 



Proof. Notice that if -My(T) < n + 1, where n G N, then 



and so, for any positive constant m, 



[ Y{s)ds <{n + l)T, 
Jo 

l^max^^s (^J^ Y{s)ds^ > m, My (T) < n + l| C {MB^Hn + l)r) > m} . (2.25) 
Using the inclusion 

{n < My{T)} C I max^Sa > (by mB), 

we obtain by ()2.25p . 

(n < My(r) < n + 1) < (^Me-Ain + l)r) > ^ - V - ^^^ ^^ ^ 

The expression for the density of the running maximum of Brownian motion \30\ Equation (2.8.4)] 
yields 



MB:,{in + 1)T) > < / ' dx. 



a 

As in [H §7.1.2], we let 



[n-y-K-dT) / {ayj {n+l)T) \/27r 



2 /"^ 2 

erfc(a) := / e""^ /^dx, Va G 
and so, 



n < My{T) < n + 1) < — erfc ' 



^ \aV(n + l)r 
Because for any a > 1, 

/■OO /-OO 

J a J a 



e-'l\ 



erfc(a) < ^e'^'/^^ Va > 1. 



we see that 



By hypothesis, T G (0,To], which implies that 

n — y — Ki9T n — y — K'&Ta 

/ > / yn G N. 

a^{n + l)T a^{n + l)ro 

Hence, provided we have 

n-y - K'dTp ^ ^ 
aV(n + l)ro - ' 
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which is true for ah n > 77,0 (y, k, cr, Tq), the smahest integer such that the preceding inequahty 
holds, we see that 



(n < My(T) < n + 1) < _^^-{n-y-n^T)y{2aHn+l)T) ^ > 



no- 



(2.26) 



Similarly, for a possibly larger no(?/, k, ??, cj, Tq), using again the fact that T £ (0, Tq], we may choose 
a positive constant c, depending also on y, k, 1?, cr and Tq, such that for all n > ng, we have 



(n — y — K'dTy 
2cj2(n + l)r ^ "2cj2r 



> c- 



Then, using the preceding inequality, we obtain the estimate (|2.24p from (j2.26p . This completes 
the proof of the claim. □ 

Next, we employ (j2.24p to obtain (|2.2ip . For any stopping time 6 G ^o,Ti we may write 



n=0 



and, by Holder's inequality, it follows 



n=0 



{A/y (T)<n+1} ^{n<MY{T)<n+l} i 



1/2 



^px(e) < ^ e^'^(^)l|MWT)<nH-i} (n < My (T) < n + 1) 



a/2 



(2.27) 



Using (|2.22p and the condition p > in (|2.20p , we have 



e l{My(T)<n+l} 
< gP(x+|r-g|T)j£2 

^ ^p(3,+|r-g|T)jg2 g2pAfsi ((n+l)T) 



exp ^2pSi y(s)dsj j l{My(T)<n+l} 

exp (^2p^max Si y(s)(is^^ l{Afy(T)<n+i} 
Vn G N (by ([iJS])). 



We see from the expression for the density of the running maximum of Brownian motion [30l 
Exercise (2.8.4)] that 



,2pA/Bi{(n+l)T) 



-x2/(2(n+l)T)^^ 



'0 ^27r(n + l)r 

< 2e2p'("+i)^, VnGN (by Mathematica), 



and so. 



-'-{A/y(T)<n+l} 



(2.28) 
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Inequalities i^i2^, (IZTTll and (l2:28]l give us 

no-l 



n=0 



(x+\r-g\T)/2 



TT 



1/4' 



in+l)T^-cn/(4.a'^T) 



no-l 



^^p{x+\r-q\T)/2 ^pHn+l)T 



n=0 



TT 



1/4 



Mx+\r-g\T)/2+p^T Y e(P^T-c/iia^T))n_ 



We choose p such that 



< p < 



2aT' 



that is, condition (|2.2U|) is obeyed, and we obtain a bound on [e^'''"^^^] which is independent of 
the choice of 6* G 5b,T- Thus, (f2:2T]l follows. (Note that ([2:21]) holds trivially when p = 0.) □ 



3. Elliptic boundary value problem 

In this section, we prove Theorem 11.31 In addition to the uniqueness result in Theorem 11.31 we 
establish the existence and uniqueness of solutions in Theorem 11.81 

The existence and uniqueness of solutions to problem (jl.ip with boundary condition (jl.3p along 
Fi, when /3 > 1, and with boundary condition (jl.Sp along dff, when < /5 < 1, are similar in 
nature. Therefore, we define 

8,<.:=^' ''"S'- (3.1) 



and treat the previous mentioned boundary value problems together as 

\ Au = f on ^, 



(3-2) 

u = g on opff. 



Now, we can give the 



Proof of Theorem \1.3\ . Our goal is to show that if u G Cioc(^U9/3^)nC^(^) is a solution to problem 
()3.2p . satisfying the pointwise growth condition ()1.20p . then it admits the stochastic representation 

We let {^fc : /c G N} denote an increasing sequence of C^"*"" subdomains of (see [25] Definition 
§6.2]) such that each tffk has compact closure in ^, and 



km 
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By applying Ito's lemma |30^ Theorems 3.3.3 & 3.3.6], we obtain for all t > 0, 



+ 1 



''y^({u,{Z{t)) + apUy{Z{t))) dWi{t) + a^l - p^Uy{Z{t))dW2{t)) . 



Since the subdomain C ^ is bounded and n G 
identity are martingales, and so we obtain 



the dVFi-terms, z = 1, 2, in the preceding 



tATff 



e-''f{Z{s))ds 



(3.3) 



We take the limit as k tends to oo in the preceding identity. By the growth estimate (j2.18p . we 
may apply the Lebesgue Dominated Convergence Theorem to show that the integral term in (j3.3|) 
converges to 

/ e-"/(Z(s))ds . 
VJo 

For the non- integral term on the left hand side of (|3.3|) , using the continuity of ?i on ^ U 5^ ^ and 
of the sample paths of the Heston process, we see that 

g-KMr^Jy^^^^ A r^J) ^ e-"(*^"*)u(Z(t A r^)), a.s. as A; ^ oo. 
Using jni Theorem 16.13], we prove that 



by showing that 



as A; —7- oo, 



is a collection of uniformly integrable random variables. By [Bl Remark related to formula (16.23)], 
it suffices to show that their p-th order moment is uniformly bounded (independent of k), for some 
p > 1. We choose p > 1 such that pMi < fi and pM2 < 1. Notice that this is possible because we 
assumed the coefficients Mi < fi and M2 < 1. Then, from the growth estimate (jl.20p . we have 

-r{tAr^,)u{z) " < Ce-'P^'^'^fc) (1 + e^^'^^ + ef*'^2^) , Vfc G N. 

From the inequality (I2.14p with c = pMi < fi and property (j2.13p applied with c = pM2 G (0, 1), 
we obtain using Mi < r /(«■!?) 

PI 



-r(tAr^,)u(^Z{t A T<?, )) < C (l + e^*'^!^ + 6^^'^^^) , Vfc G N. 



Therefore, by taking limit as k tends to oo in (]3.3p we obtain 



E^ 



-ritAre)u{Z{t A Tff))] = U{z) - Ef 



tATff 



'f{Z{s))ds 



(3.4) 



As we let t tend to oo, the integral term on the right-hand side in the preceding identity clearly 
converges to 



E^ 



'fiZ{s))ds 
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It remains to consider the left-hand side of (|3.4p . Keeping in mind that u G Cioc(^U dijff) solves 
(13.21). we rewrite this term as 



Using the growth assumption ()1.20p . we notice as above that both collections of random variables 
in the preceding identity, 

{e-'-^^<7(^(r^))lW<t} : t > 0} and {e"'-*n(Z(t))l|,^>4| : t > O} , 

are uniformly integrable, and they converge a.s. to e"^'^'^ g{Z{T^'))l^^^^^y and zero, respectively. 
Therefore, by [6l Theorem 16.13], letting t tend to oo in (|3.4p . we obtain 



□ 



[e-^^^<7(^(r^'))lw<oo}] = u{z) - / e-''f{Z{s)) ds 
which implies that u = u* on & \Jd^G ^ where u* is defined by (jl.23p . 



Proof of Theorem\r^ Our goal is to show that if < /3 < 1 and it e Cioc(^ U Ti) n C2(^) n 
C^'[^^(^U To) is a solution to problem p.ip . satisfying the growth estimate (|1.20p . then it admits 
the stochastic representation (ll.25p . 

We consider the following sequence of increasing subdomains of ff, 



'^/k := {z £ : \z\ < /c,dist (z^Ti) > 1/k} , /c G N, 
with non-empty boundary portions Tq D '^k- Let e > and denote 

:=¥ + e, and := (X, Y^) . 
By applying Ito's lemma [30^ Theorems 3.3.3 &: 3.3.6], we obtain 



E? 



e-''A'u{Z'{s))ds 



yt > 0, 



where z/o^^. is given by ()1.22p . and denotes the elliptic differential operator, 

S S r 

A^v := Av + -Vx + KSVy - - {va:x + ^pav^y + cf'^Vyy) , yv e 
Using (II. ip . we can write (13. 7p as 



(3.5) 
(3.6) 

(3.7) 
(3.8) 



E^ 



U[Z 



E? 



e-''fiZ%s))ds 



El 







e-''%A'' - A)u{Z'{s))ds 



(3.9) 



First, we take limit as e tends to in the preceding identity. We may assume without loss of 
generality that e < 1/k, for any fixed A; > 1. We evaluate the residual term {A^ — A)u with (j3.8p 
to give 

(3.10) 



\{A'-A)u{Z%s))\<Ce\\Du\\c^^^^^+C(^l^y^^^^^^^^ + ^j \\yD'u\\c 

for all < s < t A I'ap/,,, where C is a positive constant depending only on the Heston constant 
coefficients. This follows from the fact that 



eDMZ%s))=eD\{Z%s))l^y,^^^^^^^+eDMZ%s))l^y^^^^^^^y 



Vs > 0, 
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Combining the preceding inequality with the definition (j3.8p of A^, we obtain (|3.10p . Since u G 



C. 



s.loc ^ 



' U To), and 



l{y^'(s)<v/i} ^ 0' aseiO, 
we see that by p.lOp yields, for each k > 1, 



e-''{A' - A)u{Z'{s))ds 



0, as e I 0. 



(3.11) 



In addition, using the continuity of / and u on compact subsets of ^ U Fq, we have 



e-'V{Z'{s))ds 







e-'''f{Z{s))ds 





as e I 0, 
as e I 0. 



(3.12) 



Therefore, using (j3.1ip and the preceding limits, we find that p.9p gives 



-(^^■^^,)u{Z{tAu,,J)] =u{z)-Ef 



'f{Z{s))ds 



Note that by letting k and t tend to oo, we have 

t A z^f/j. — )■ i^^^, a.s. 



(3.13) 



(3.14) 



By using the same argument as that used in the proof of Theorem 11.31 to take the limit as k and t 
tend to oo in (|3.3p . we can take the limit as k and t tend to oo in (|3.13p to give 



[e-''"^g{ZM)]=u{z)-EI^ 



^0 



e-''^f{Z{s))ds 



This establishes u = u* , where u* is given by (|1.25p . and completes the proof. 



□ 



Next, we prove existence of solutions to problem (j3.2p when the boundary data g is continuous 
on suitable portions of the boundary of ^. 



Proof of Theorem Following the comments preceding problem (|3.2p , we need to show that u* , 
given by (fT^ . is a solution to problem ([22]), that u* e Qod^ U dpff) nC'^^ff), and that u* 
satisfies the growth assumption (|1.2Up . 

Notice that Lemma |2.12| applied with 6i = t^, 02 = oo and -0 = 0, shows that u* defined by 
fjl.23p satisfies the growth assumption (jl.20p . It remains to prove that u* £ Cioc(i^U5/3^)nC^(^). 
Notice that Theorem 11.31 implies that u* is the unique solution to the elliptic boundary value 
problem (13. 2p . since any C\oc{^ U dj^ff) (1 C'^{ff) solution must coincide with u*. 

By hypothesis and the definition of dsff in (|3.ip . we have g S Cioc( 



). Since is closed, we 
may use j21|, Theorem 3.1.2] to extend g to such that its extension g G Cioc(M^). We organize 
the proof in two steps. 
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Step 1 {u* G C^+°((^)). Let {^k ■ k e N} be an increasing sequence of C^+" subdomains of 
as in the proof of Theorem II. 3[ We notice that on each domain the differential operator A 
is uniformly elliptic with C°°(ffk) coefficients. From our hypotheses, we have / G C^{i^k) and 
g G C(^fc). Therefore, Theorem 6.13] implies that the elliptic boundary value problem 



(3.15) 



admits a unique solution u^. G 
stochastic representation on &k 




Moreover, b}0 [Ml Theorem 6.5.1], Uk admits a 



Uk{z) = E| 



'=5(^(-^<>fc))l{r^^<oo} 





[/ 




Jo 



'f{Z{s))ds 



(3.16) 



Our goal is to show that converges pointwise to u* on ^. Recall that r^, is an increasing 
sequence of stopping times which converges to almost surely. Using g G C\oc{^ U dpi^) and 
the continuity of the sample paths of the Heston process, the growth estimate (jl.20p and Lemma 
\2.12\ the same argument used in the proof of Theorem 11.31 shows that the sequence {uk ■ k G N} 
converges pointwise to u* on 0'. 

Fix zq := (xo,2/o) G and choose a Euclidean ball B := B{zQ,ro) such that B d ff. We denote 
i?i/2 = B{zQ,ro/2). As in the proof of Lemma [2.12^ the sequence Uk is uniformly bounded on B 
because it obeys 



yz = {x,y) eB,ke N. 



From the interior Schauder estimates |25t Corollary 6.3], the sequence {uk ■ k G N} has uniformly 
bounded C^''~°(-Bi/2) norms. Compactness of the embedding C^^"(-Bi/2) C'^'^'^ for 
< 7 < a, shows that, after passing to a subsequence, the sequence {uk : /c G N} converges in 
C'^^'^ {B112) to u* G C'^^'^ {Bii2)-, and so Au* = f on Bi/2- Because the subsequence has uni- 
formly bounded C'^~^'^{Bi/2) norms and it converges strongly in C'^{Bi/2) to n*, we obtain that 

U* G C2+"(Si/2). 

, SO it remains to 



Step 2 {u* G Cioc(i^U dpi?)). From the previous step, we know that u* G 
show continuity of u* up to dpff. We consider two cases. 

Case 1 (n* G Cioc(^U Fi), for all (3 > 0). First, we show that u* is continuous up to Fi. We fix 
2o G Fi, and let B be an open ball centered at zq, such that B n 9IHI = 0. Let U := Bf^ff. Let the 
function g be defined on dU such that it coincides with g on dU H Off, and it coincides with u* on 

^ur^&. 



Claim 3.1. The strong Markov property of the Heston process {Z{s))s>o CLnd the definition ()1.23p 
of u* , implies that 



u*{z)=E^Q [e---^g{Z{Tu))]+El 



-rt 



f{Z{t))dt 



\/z G U. 



(3.17) 



Proof. By Corollary 12. 8| the Heston stochastic differential equation (jl.l7p admits a unique strong 
solution, for any initial point {t,x,y) G [0, 00) x M x [0, 00), and [16l Theorem 1.16(c)] shows that 
the solution satisfies the strong Markov property. 



^See also [SO] Proposition 5.7.2], [38l Theorem 9.1.1 & Corollary 9.1.2]. 
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Let z £ U, then T^ < a.s. Since Z is a time-homogeneous strong Markov process, we obtain 



-rrxj: 



which can be written as 



TU='l'ff}\ 



+ E? 



{Ti7<T^} 



(3.18) 



Similarly, we have for the integral term 



E^ 



re 



-rt 



f{Z{t))dt 





\l 




Jo 



-rt 



f{Z{t))dt 



■El 



^{tu<T0} 



-rt 



f{Z{t))dt 



and so, by conditioning the second term in the preceding equality on ^(tjj) and applying the 
strong Markov property, we have 



El 



'-{ru<r0} 



-rt 



f{Z{t))dt 



El 



El 



Hru<rff} I e 
'ru 



f 



-rt 



{ru<TfyY 



f{Z{t))dt 

T0 



(3.19) 



Combining (l3J8]l and (l3J9]l in ([L23]), we obtain 

= E^ [e~^^^'giZiTu))l{ru=re}] + n 



-rt 



f{Z{t))dt 



+ EI. 



T0 



-rt 



f{Z{t))dt 



Using again (I1.23P for u*{Z{tu)), the preceding equality yields (I3.17p . This completes the proof of 
Claim EH □ 

By j251 Theorem 6.13] and a straightforward extension o^ \2A\ Theorem 6.5.1] from domains 
with to domains with regular boundary, as in |12t §6. 2. 6. A], the integral term in ()3.17p is the 
solution on U of the uniformly elliptic partial differential equation Au* = f with homogeneous 
Dirichlet boundary condition, and it is a continuous function up to dU. Notice that dU satisfies 
the exterior sphere condition and thus dU is regular bj0 [la Theorem 2.4.4] (see \12\ Definitions 
2.4.1 & 6.2.3] for the definitiorH of regular points of dU). The continuity of the non-integral term 
in ()3.17p at zq follows from Corollarv IB. 2) as g is continuous at zq by hypotheses. 

It remains to show that, when < /3 < 1, the solution u* is continuous up to Tq. 

Case 2 [u* G Cioc('^U fo), for all < /3 < 1). Let zo = (xo,0) E fo. We denote by 6*^ the first 
time the process started at z = {x, y) £ & hits y = 0. Obviously, we have 

Tff<e'< a.s., (3.20) 

where Tq is given by (j2.1ip . For /3 G (0, 1), it follows from (j2.12p and the preceding inequality 
between stopping times, that 9^ converges to 0, as y goes to 0, uniformly in x S M. Therefore, 



^See also [Ml Proposition 5.7.2], [38l Theorem 9.1.1 & Corollary 9.1.2]. 
^See also [Ml Proposition 4.2.15 & Theorem 4.2.19]. 
5c 



^See also [30l Definition 4.2.9], [38l Definition 9.2.8]. 
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(3.21) 



the integral term in (|3.17p converges to zero. Next, we want to show that the non-integral term in 
(|3.17p converges to g{zo). We rewrite that term as 

+ <7(^o)(l-IE^ [e"'"''])- 
From the observation that r|. < 9^ a.s., we see that 

[e"''^^] ^1, as z ^ zq. (3.22) 

By (j3.2ip . it remains to show that '&^[e~''"^^ {g{Z{Tff)) — q^zq))] converges to zero, as z E ^ 
converges to zq. We fix e > and choose 5i > such that 

\g{z)-g{zo)\<e, Vz G S(zo, (^i) n 9^. (3.23) 

From [301 Equation (5.3.18) in Problem 5.3.15 ], there is a positive constant Ci, depending on zq 
and 5i, such that 



sup E^ 



sup \Z{s) — z\ 

0<s<t 



from where it follows 



sup 

z£B{zo,5i 



sup \Z{s) -z\> 6i/2 < 

0<s<t ) 0\ 



Next, we choose t > sufficiently small such that 

2CWt 

and, by (j3.20p and (j2.12p . we may choose bi> ^ sufficiently small such that 



< e, 



Tf:' > t ) < e. 



(3.24) 

(3.25) 
(3.26) 



Let (5 := min{5i/2, (52}- We rewrite 



[g[Z(Te)) - g[z^)) = e-'-^^ (^(^(r^)) - 5(^o)) \r„<t) 

{g{Z{Tff)) - g{zo)) l{r^>t} 



+ e 



to give 



{g{Z (Tfj)) - g{zo)) = e '"^^ {g{Z (rff)) - g{zo)) l{r^<f ,supo<,<, \z(s)-z\<Si/2} 



+ e {g{Z(Te)) - gizo)) l{r^<Mupo<,<i \z{s)^z\>5i/2} 
+ e-'^^ ig{Z{Tff))-gizo))l{r^:.t} 



By (I323I), we have for all z £ B{zo, 5) n 



E^ 



\giZ{Tff)) - g{zo)\ l{rff<t, supo<,<i \Z{s)^z\<5i/2} 



< e. 



(3.27) 



(3.28) 



We choose p > 1 such that pA-Ii < fi and pM2 < 1. Notice that this is possible because we assumed 
the coefficients Mi < ^ and M2 < 1. Then, from the growth estimate p.20p for g, we have 

From the inequality (|2.14p with c = pMi < fi and property (j2.13p applied with c = pM2 £ (0, 1), 
we obtain using the condition Mi < r/{K'd) 

[\e-'^^g{Z{Tg))\''] <C{1 + eP^^y + e^'^^^) . 
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Let C2 > be an bound on the right-hand side of the preceding inequahty, for ah z = {x, y) G 
B{zq,6) n (y. By the Cauchy-Schwarz inequahty, we have 

[e— ^ {g{Z{Tg)) - gizo)) l^r^yt}] \ < [e"''^^'' |5(^(t^0) " ^(^o)^] '^'^ (r^^ > tf"' , 

where p' > 1 denotes the conjugate exponent of p. Using the fact that t|. < Tq^ from ()3.20p and 
()3.26p . we obtain in the preceding inequahty 

[e— {g{Z{T,^)) - g{zo)) l{.,>t|] | < 2C',/'Q^ (To > tf^' 

< 2Cl^^e^/P', Vz G B{zo, 6) n ff, 

From the inequahty, 



(3.29) 



{g{Z{Tff)) - 5(2:0)) l{r«?<i,suPo<.<t \Z{s)-z\>Si/2} 



1/p' 



< [e-^P-^ \g{Z{T^))-g{zor]'/''Q' ( sup |Z(s) - z\ > 6^/2 

\0<s<t 

the inequahties ()3.24p and ()3.25p and definition of C2 yield 

[e^''"'' i9{Z{T^)) - g{zo)) lw<Mupo<.<, \zis)-z\>s,/2}] I < 2C2'/V/f'. (3.30) 
Substituting (lOHD . (lOOD . and (lOOD in ([3211), we obtain 

E [e-^-^ {g{Z{Tff)) - 5(zo))] < (1 + ^Cl'") (e + e^/^') , G i?(zo, 5) n 
and so n* is continuous up to Fq, when < /3 < 1. 

This concludes the proof that u* £ C\oc{G yjdpff), for ah ^5 > 0. 

This completes the proof of Theorem 11.81 □ 

We now prove existence of solutions to problem ()3.2p when the boundary data g is Holder 
continuous on suitable portions of the boundary of i?. 

Proof of Theorem \l.y[. The proof of the theorem is the same as that of Theorem 11.81 with the 
exception that Case [H of Step [2] can be simplified by applying the classical boundary Schauder 
estimates. Also, instead of using the sequence of subdomains {i^fc : A; G N} precompactly contained 
in ff, as in the proof of Theorem ll.3| we consider an increasing sequence, : A: G N}, of C^"*"" 
subdomains of (see [25l Definition §6.2]) such that each &k satisfies 

(-A:,A:) X (l//c,A:) C ^fc C (-2A;,2/c) X (l/(2/c),2A:), V/c G N, (3.31) 

and 

U = ^• 

fceN 

Since Fi is assumed to be of class C^"*"", we may choose Sik to be of class C^+". 
Let £ Ti, and tq > small enough such that B{zQ,rQ) n Fq = 0. Let 

D := B{zo,ro) D G and D' := B{zo,ro/2) D G. 

By (I3.3ip . we may choose /cq G N large enough such that D C ^fc, for all k > /cq. Using / G C'^{D), 
g G C^"^"(Z)) and applying [25] Corollary 6.7], and the fact that Uk solves (|3.15p with g replacing 
g and replacing G^, we have 

Il^^fc|lc2+"(B') ^ C' (^||'Ufc||c{D) + ll5'llc2+"{5) + ll/llc"(S)) 1 ^ ^0, (3.32) 
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where C > is a constant depending only on the coefficients of A, and the domains D and 
D' . Combining the preceding inequaUty with the uniform bound on the C{D) norms of the se- 
quence {ufc : k G N}, resulting from Lemma [2 .121 the compactness of the embedding of C'^^"^{D') 
C'^^^[D'), when < 7 < a, implies that a subsequence converges strongly to u* in C'^~^^{D'). 
Therefore, u* G Cl+^{D'), and Au* = f on D' and u* = g on dD' D Ti. Moreover, u* G C^+^fl)'), 
since Uk G C'^~^"'{D'), for all k > ko, and the sequence converges in C^+'fiD') to u*. Combining the 
boundary estimate (I3.32p with Step [U in the proof of Theorem 11.81 we obtain the conclusion that 

n* G cf+"(^uri). □ 

Remark 3.2 (Validity of the stochastic representation for strong solutions). The stochastic rep- 
resentation (|1.25p for solutions to problem (jl.ip with boundary condition along Fi is valid if we 
replace the condition u G Ciod^ U Ti) n C^{ff) D Clf^^i^ U Tq) in the hypotheses of Theorem (LSI 
with the weaker condition u G Ciod^ U Ti) n W^^^{^) H Cl'l^^ U Tq). 

4. Elliptic obstacle problem 

This section contains the proofs of Theorems 11.121 and 11.131 As in problem (j3.2p , the questions of 
uniqueness of solutions to problem (jl.2p with Dirichlet boundary condition along Fi, when /3 > 1, 
and along dff, when < /3 < 1, are similar in nature. We can conveniently treat them together as 



min {Au — f,u 
u = g 



on 
on d. 



'pi 



(4.1) 



where dpff is given by ()3.ip . 

Proof of Theorem \1.12\ Lemma [2. 121 indicates that u* given by ()1.27p satisfies (11.20p . so the growth 
assumption on u in Theorem 11.121 is justified. 

By the preceding remarks, it suffices to prove that the stochastic representation ()1.27p holds for 
solutions u G Cioc(^U 5/3^) fl C^(^) to problem ()4.ip . We consider the two situations: u > u* and 
u<u* on ffUdi30', where u* is defined by (fL271) . 

Step 1 (Proof that u > u* on & U dpi?). Let {i^a: : A: G N} be an increasing sequence of C^"*"" 
subdomains of as in the proof of Theorem 11.31 Since u G C^(^), Ito's lemma j30| Theorems 
3.3.3 & 3.3.6] yields, for any stopping time G =^5^, 



MeAr^,)u^Z{e A T^,))l = uiz) - E, 



6»At 



e-'''Au{Z{s))ds 



(4.2) 



By splitting the right-hand side in the preceding identity. 



-(^^-^.)t.(Z(eAr^g)' 



and using u> ip on and Au > f a.e. on i^, the identity 



■%(Z(0Ar^J)l|, 
T2]) gives 



uiz) > Ef 



-re 



+ Ef 



'f{Z{s))ds 



(4.3) 



Just as in the proof of Theorem II. 3|, the collections of random variables 



{e-^-U{z{e))i 



: A: G N 



I and |f 



: A: G N 
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are uniformly integrable because u and ip satisfy the pointwise growth estimate (jl.20p . From the 
continuity of u and ip on Or^, we also have the a.s. convergence, 



-re 



} 



-re 



ip{Z{e))l{0^^^y, as k 



oo, 



e ""*fcn(Z(r^.J)l|,^^<e} ^ e u(Z(r^))l|,^<e}, as /c ^ oo. 

Therefore, by [51 Theorem 16.13], we can take limit as k tends to oo in inequality ()4.3p to obtain, 
for all 9£ 3^, 



■Eh I e-''fiZis))ds 



which yields u> u* on L) djjff. 

Step 2 (Proof that u < u* on U dpi)"). The continuation region, 

<r := {u > V'}, (4.4) 

is an open set by the continuity of u and V- We denote the first exit time of Z*'^ from the 
continuation region, by 

f*'^ := {s > t : Z*'^(s) ^ ^} , (4.5) 

and write f = f*'^ for brevity. This is indeed a stopping time because the process Z*'^ is continuous 
and ^ is open. By the same argument used in Step [1] with 6 replaced by f, we obtain that all 
inequalities hold with equalities because u{Z{t)) = ip{Z{f)) and Au = / on the continuation region. 
Therefore, 



(z) = [e-^^^(Z(f ))1|,<,^}] + E^ [e-^^-^<7(^(r^0)lw<r>}] 
e~^'f{Z{s))ds , 



+ EI 



which implies that u < u* . 

By combining the preceding two steps, we obtain the stochastic representation ()1.27p of solutions 
to problem ()4.ip . and hence the uniqueness assertion. □ 

Proof of Theorem \1.13[ Lemma [2. 121 indicates that u* given by (jl.28p satisfies (jl.20p . so the growth 
assumption on u in Theorem 11.121 is justified. 



Our goal is to show that if < /3 < 1 and u G Cioc{^LlTi)nC^{0')nC, 



1,1 ( 



'"UTo) is a solution to 

problem (|1.2p with Dirichlet boundary condition (|1.4p along Fi, and satisfying the growth estimate 
(|1.20p , then it admits the stochastic representation (|1.28p . As in the proof of Theorem I1.12| we 
consider the following two cases. 

Step 1 (Proof that n > n* on U Pi). Let e > and : k G N} be the collection of increasing 
subdomains as in ()3.5p . By applying Ito's lemma, we obtain, for all t > and 9 € 



u{z) = E'q re^"(*^'^*fe^^)n(Z"(t A i/.^, A 6)) 



'A'uiZ'{s))ds 



(4.6) 
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where 1^94 is given by (|1.22p and is defined in (j3.6p . and is defined by (|3.8p . By (jl.2p and 
(|3.8p . preceding identity gives 



n(z) > [e~"(*^^*fe^^)n(Z^(t A z^.^, A 9)) 



'f{Z'{s))ds 





[L 







e-''\A' - A)u{Z'{s))ds 



(4.7) 



First, we take the limit as e tends to in (j4.7p . We can assume without loss of generality that 
£ < 1/k, for any fixed k € N. The residual term (A^ — A)u then obeys estimate (jS.lOp because 
''UFo). Therefore, (|3.1ip also holds in the present case. In addition, using the continuity 



s,loc^ 



of /, u, Du and yD^u on compact subsets of ^ U Tq, we see that (13.12P holds, and so, by taking 
hmit as e I in (fiTp . 



u{z) > El 



(4.8) 



Finally, letting k and t tend to 00 and using the convergence (j3.14p . the same argument employed 
in the proof of Theorem 11.31 can be applied to conclude that u > u* on 0' UTi, where u* is given 
by (fL2H]l . 

Step 2 (Proof that u < u* on ^ U Fi). We choose 9 = f in the preceding step, where f is defined 
by ()4.5p . By the definition ()4.4p of the continuation region, 'rf, and the obstacle problem ()1.2p . we 
notice that inequalities ()4.7p and (|4.8p hold with equality and so it follows as in Step [J that u < u* 
on i^U Fi. 



This completes the proof. 



□ 



Remark 4.1 (Validity of the stochastic representation for strong solutions). The stochastic repre- 
sentation (|1.27p of solutions to problem (j4.ip . when /3 > 0, holds under the weaker assumption 



that u G C] 



loci 



U5« 



2,2 

loc 



Similarly, the stochastic representation (jl.28p of solutions to 



problem (jl.2p with Dirichlet boundary condition ()1.4p along Fi, when < /3 < 1, holds under the 



weaker assumption that u G Cioc(i^UFi)nC, 



s,loc^ 



UFo)UH^, 



2,2 
loc 



In each case, we would replace 



the application of the classical ltd lemma |30( Theorems 3.3.3 & 3.3.6] with [5l Identity (8.62) in 
Theorem 2.8.5], or we could apply an approximation argument involving C^(^) functions. 



5. Parabolic terminal/boundary value problem 

This section contains the proofs of Theorems 11.161 and 11.181 and an existence result in Theorem 
11.221 Because the Heston process satisfies the strong Markov property, it suffices to prove the 
stochastic representation of solutions to the terminal value problem for T as small as we like. In 
particular, without loss of generality, we can choose T such that 

Hypothesis 5.1. There is a constant po > 1 such that 

(1) Condition (|2.20p in Lemma [2.14l is satisfied for p := pqM2, where M2 G [0, 1] is the constant 
appearing in (jl.29p : 

(2) One has poMi < fi, where Mi G [0,^) in (fOal) . 

As in ^ we first prove uniqueness of solutions to the parabolic terminal/boundary value prob- 
lems (|1.7p with different possible Dirichlet boundary conditions depending on the parameter /3. 
The proofs are similar those of Theorems 11.31 and 11.51 
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The existence and uniqueness of solutions to problem (jl.7p with boundary condition (jl.9|) . when 
/? > 1, and with boundary condition (jl.l2p . when < /3 < 1, are similar in nature. By analogy 
with our treatment of problem (j3.2p . we define 



(5.1) 



gig if/3>i, 

SQ ifO</3<l, 

where we recall that Q := (0, T) x ^. The preceding problems can then be formulated as 

-ut + Au = f on Q, (5.2) 
u = g on <5pQ. (5-3) 

We now have the 

Proof of Theorem \1.1(A We choose T > small enough and pq > 1 as in Hypothesis 15.11 The 
pattern of the proof is the same as that of Theorem 11.31 For completeness, we outline the main 
steps of the argument. 

We need to show that if u G Cioc('3 U S/3Q) ^ C^iQ) is a solution to problem (j5.2p . satisfying the 
growth bound (jl.29p . then it admits the stochastic representation (|1.30p . We choose a collection 
of increasing subdomains, ■ k S N}, as in the proof of Theorem 11.31 By applying Ito's lemma 
|30^ Theorems 3.3.3 & 3.3.6], we obtain, for all t > and G N, 



t,z 



u{t,z) - E, 



t,z 



-'■^'-'^f{s,Z{s))ds 



(5.4) 



We now take limit as k tends to 00 in the preceding identity. Using (jl.29p and Lemma 12.111 we 
obtain 



E, 



i,2 



-^(^-'^f{s,Z{s))ds 



E 



t,z 



-^('-'^f{s,Z{s))ds 



as k 



00. 



(5.5) 



From the continuity of u and of the sample paths of Z, we obtain the a.s. convergence as k tends 
to 00, 

^-rir^^AT-t)^^^^^ A T, Z(t^, A T)) ^ e^"-^-^^^^ g{r^^ A T, Z(t^ A T)). 
In order to prove that, as k tends to 00, 



e: 



t,z 



-<^^k^^-')u{T^^AT,ZiT^^AT)) 



E, 



t,z 



Mre^T)g^^^,^r^^ZiT^AT)) 



(5.6) 



using [6l Theorem 16.13], it is enough to show that the collection of random variables, 

|e-r{r^,AT-t)^(^^^ A T, Z(t^, A T)) : A; G n} (5.7) 

is uniformly integrable. For po > 1 as in Hypothesis 15. 11 it is enough to show that their po-th. order 
moments are uniformly bounded ([U Observation following Equation (16.23)]), that is 



sup E, 

fceN 



t,z 



-rre. 



(r<y^,Z(T^J)l| <T} 



po 



< +00. 



(5.8) 



From (jl.29p . we have, for some constant C, 



E, 



t,z 



-'i-'^k^^-'^u{Tff^AT,Z{T^^AT)) 



< C 1 + E, 



t,z 



.PoMjYiTff.AT) 



Po 



t,z 



^PoAhXiTff AT) 
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Now, the uniform bound in (|5.8p follows by applying the supermartingale property (j2.14p with 
c := pqMi to the first expectation in the preceding inequality, and by applying (j2.2ip with p := P0M2 
to the second expectation above. Therefore, by taking the limit as k tends to 00 in (j5.4p . with 
the aid of (I5.5p and (15 .Gp . we obtain the stochastic representation ()1.30p of solutions to problem 
(|^ . □ 



Proof of TheoremlTM The need is to show that if < /3 < 1 and n G CiodQ U d^Q) H C^{Q) D 
Cg'[]^^((0, T) X (i^UFo)) is a solution to problem ()1.7p with boundary condition (|1.9p . satisfying the 
growth bound (ll.29p . then it admits the stochastic representation (ll.32p . 

Let e > and : A; € N} be the collection of increasing subdomains as in (13.5p . By applying 
Ito's lemma, we obtain 



E, 



t,z 



~-(TAu^^Ju{T A i/.^, , Z%T A z..^J)l = u{t, z) - E, 



e-'^'A^u{s,Z^{s))ds 



where vo^^ is given by (jl.22p . Z'' by ()3.6p and is defined by ()3.8p . The proof now follows the 
same path as that of Theorem 11.51 with the only modification being that we now take the limit as 
k tends to 00 in the preceding identity in order to obtain (jl.32p . □ 



Analogous to Lemma 12.121 we have the following auxiliary result. 

Lemma 5.2. Suppose f and g obey the growth assumption (11.29p . Then there are positive constants 
C , Ml < (JL and M2 € [0, 1], such that for any stopping times 61,62 G c%^t with values in [t, T], the 
function Jp^'^^ given by (I1.36P obeys the growth assumption (ll.29p . 

Proof. The proof follows as in Lemma 12.121 with the aid of Lemma 12.111 Notice that because the 
stopping times 61,62 G =3^,t are bounded by T, we do not need the constant r to be positive, as in 
Lemma [2I2I ' □ 

Remark 5.3. The function ^ in ()1.36p plays the role of the obstacle function and is relevant only 
for problem (jl.Sp . 

Next, we have the following existence results for solutions to the parabolic terminal/boundary 
value problem ()5.2p . for all /3 > 0. 



Proof of Theorem \1.22\ We choose T > small enough and pq> 1 as in Hypothesis 15.11 

By hypothesis, we have g G Cioc(S/3Q)- Since QpQ is closed, we may use [211 Theorem 3.1.2] to 
extend 5 to a function on [0,T] x M^, again called g, such that g G Cioc([0,T] x M?). 

The proof follows the same pattern as that of Theorem 11.81 For completeness, we outline the 
main steps. Let be an increasing sequence of C^"*"" subdomains of & as in the proof of Theorem 
11.31 and let Qj. := (0,T) x 0'^. We notice that on each cylindrical domain, Qk, the operator A is 
uniformly elliptic, and its coefficients are C°°{Qk) functions. By hypothesis, there is an a G (0, 1) 
such that / G C°(Qfc) and g G C{Qk)- Therefore, by [211 Theorem 3.4.9], the terminal value 
problem 

-Ut + Au = f on Qk, 

u = g on (0,r) X 9^fcU{r} X ^fc. 
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has a unique solution Uk G C{Qk) H C^+"((5fc), and [Ml Theorem 6.5.2] it has the stochastic 
representation 

'■rff^^^T 



-<'-'^f{s,Z{s))ds 



(5.9) 

V(t,z) G Qfc. 



Because r^^. converges to a.s. as /c — )■ oo, the integral term in (j5.9p converges to the integral 
term of u* in (jl.30p , by the same argument as that used in the proof of Theorem 11.161 By the 
continuity of g and of the paths of the Heston process Z, we also know that 

^-r(r„^AT)g^^^,^ A T, Z{Tff^ A T)) ^ e-^^^^^^)^!^^' A T, Z{Te A T)), as ^ oo. 

In order to show that the preceding convergence takes place in expectation also, it is enough to 
show that the collection of random variables, 



|e--(-^.^^)ff(r^^ A T, Z^r^^ A T)) : A; G n} 



is uniformly integrable, but this follows by the same argument as that used for the collections 
(j5.7p in the proof of Theorem 11.161 by bounding their pQ-lh. order moments (po > !)• Therefore, 
the sequence {uj. : k G N} converges to u* pointwise on Q. By interior Schauder estimates for 
parabolic equations [24^ Theorem 3.3.5] and Lemma 15.21 there is a subsequence of {uk ■ k G N} 
which converges to u* in C^''~"'((5), when < a' < a. Using the Arzela-Ascoli Theorem, we obtain 
u* G C'^~^'^{Q). The proof of continuity of u up to S^Q follows by exactly the same argument as 
that used in the proof of Step [2] in Theorem 11.81 Therefore, u* is a solution to ()5.2p . 

From Theorem 11.161 and Lemma 15. 2| we see that u* in (ll.30p is the unique solution to the 
parabolic terminal value problem ()5.2p . for all /3 > 0. □ 

We now have the 

Proof of Theorem \1.23[ Just as in the proof of Theorem 11.221 we can easily adapt the proof of 
Theorem 11.91 for the elliptic case to the present parabolic case. For this purpose, we only need to 
make use of the local boundary Schauder estimate in Proposition lA. II instead of [251 Corollary 6.7] 
for the elliptic case. □ 

6. Parabolic obstacle problem 

Problem (jl.Sp with boundary condition (jl.l2p . when < /3 < 1, and with boundary condition 
(II. 9p . when /3 > 1, can be formulated as 



I min {-tij + ^ti - /,n -■(/'} = on Q, ^^^^ 
\u = g on g^Q, 

where S/jQ is defined in (|5.ip . According to Theorem 1 1.27| the solution to problem (|6.ip is given 
in (fL38D . 



Proof of Theorem \1.27\ We choose T > small enough so that it obeys Hypothesis 15.11 For such 



T > 0, the proof of Theorem 11.121 adapts to the present case in the same way that the proof 
of Theorem 11.31 adapts to give a proof of Theorem 11.161 Therefore, it remains to show that the 
corresponding stochastic representation (|1.37p of the solution to problem (16. ip holds for T arbitrarily 
large. 



^See also [SO] Theorem 5.7.6]. 
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Let N := [T/T\ (the greatest integer in T/T), and Tj := iT, for i = 0, . . . , iV — 1, and T/v := T. 
Let k be an integer such that 1 < k < N — 1, and assume that the stochastic representation formula 
(|1.37|) holds for any t G [TJ, T], where i = k, . . . , N — 1. We want to show that it holds also for 
t G [Tk-i,T]. Notice that for k = N — 1, we have T — t < T, for all t G [Tjv„i, T], and so we know 
that the stochastic representation (I1.37P of the solution to problem (ll.Sp holds, by the observation 
at the beginning of the present proof. 

For any t < v < T, stopping time 6 G c%^v with values in [t,v], and G C(^), we denote 

Jt [0.2) 
Notice that by choosing = g{T, •) and u = T in (j6.2p . we obtain, for any stopping time 9 G c%^t, 

and so, 

+ e-^-(^^^^-*)<7(r^^ A T, Z(t^ A T))l|,^;,r<e}. 

Because u solves problem (|6.ip on the interval (Tfc_i,Tfc), and — Tfc_i < T, we see that u has 
the stochastic representation (ll.37p . for any t G [Tfc^ijTfc) and z G 9/3^, 

n(t,z)= sup eJ;^ [F"*(^'--)(i,z,r,.,0)l . (6.4) 
For any stopping time r\ G =3^,Tfe > we set 

e-'^'-'^f{s,Z{s))ds 
+ e-^(^-*)^(r?,Z(r?))l|,<,^;,^„l (6-^) 

+ e-^("'^-*)5(T^, ^(T^O)lW<r,,r,<T,.}, 

and for any stopping time ^ G 3^Tk,T, we let 

+ e-^«-^'=)V'(e,^(e))i{c<r,AT} 

+ e-^(^'^^^-^^)5(r^^ A T, Z(t^ A T))1{^^^t<0- 
For the rest of the proof, we fix z G & U dpff and i G [Tfc_i, T^). 
Let rj G =!>^,Tfc and ^ G '^Tk,T- It is straightforward to see that 

\C if ^ = Tk, 
is a stopping time with values in [t,T]. We denote by 

yt,T = G ■^t,T ■■ = ^1{„<T,} + a{„=Tfe}' where r/ G 5J,t, and ^ G ^t„t} ■ (6.7) 



(6.6) 
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For any stopping time 9 G =3^,Ti we define the stopping times 6' G ^t,Tu 9" G ■^Tk,T, 

^' ■= l{6»<Tfc}^ + ^{e>n}Tk and 9" := l{e<j.^^}Tfc + l{6»>Tfc}^'- 
Then, any stopping time 9 G c%^t can be written as 

9 = 9'lf^Q^Tk} + ^"^{e>n} 
= 6''l{e'<Tfc} + 6'"l{e'=Tfe} 

and so, 

The preceding identity and definitions ()1.37p of u* and ()6.2p of F''' give us 

w*(t,z)= sup E^;^ [F^^(^'-)(t,z,r,0) 

We shall need the following identities 



(6.9) 



Claim 6.1. For any stopping time 9 = 'i]l^^^Xk} ^'^{v=Tk}' ^^^re i] G '^t,Tt 0.''^^ ^ G 3^Tk,T, we 
have the following identities 

Tghd fTfff\r) 

e-^'(^-*)/(., Z{s))ds = l|,<Tj j e-'^'-'^f{s, Z{s))ds 



+ 1{,=T,} r'^^ e-'-^'-'^f{s,Z{s))ds, 



and 



and 



-<'-'H{9,Z{9))l{e<rg;,T} = e~^^^~'^{r^,Z{r^))l{^^,^^r,}l{r,<n} 

+ e-^«-*)V(€,^(C))l{e<r^AT}l{,=r,}, 



Proof. Notice that 

{9 < n} = {rj< Tk} and {9 > = {rj = TJ. (6.10) 
The first identity is obvious because, by (|6.10p . we see that 

9 = r] on {rj < T^} and 9 = on {r] = T^}. (6.11) 
The second identity follows by the observation that 

{9 < Tff ^T} = {9< re NT,9 < Tu} U{9 <t^AT,9> U}, 
and using (16. lip and (|6.10p . it follows 

{9 < Tff A T} = {?? <tgN Tk,v < Tk} U < T^^ AT,7] = Tk}. 
For the last identity of the claim, we notice 

{t^ AT<9} = {TffAT<9, Tff <T}U {Tff AT<9,Tff> T} 

= {t0 at <9,Tff <T,9 < Tk} U {t^ at <9,tg <T,9> Tk} 
U{Tff AT <9,TG>T}. 
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By ([6TTD and IKTOh . we obtain 

{T^AT<e} = {Te <r],T0<T,T]< Tk} U {t^ A T < ^,Tff <T,t] = Tk} 
U{r^^AT<e,r^' >T} 
= {Tff<r],r]< Tfc} U{TffAT<^,r] = U}, 
which imphes the last identity of the claim. □ 
We can write the expression for F3^'^''\t, z, T, 9) as a sum, 

F5(^'-)(t,z,T,0) = l{,,^T,}F^{t,z,Tu,ii) + l{.,=T,}e-'^'^'~'^F2{t,z,T,i). (6.12) 

Because ^ G ^Tk,T and F2{t, z,T,^) depends only on (Z*'^(s))^^^^^^, and the Heston process has 
the (strong) Markov property Theorem 1.15 (c)], we have a.s. that 

E'A^[F2{t,z,T,CWn]=Ef^"^'-^-^ [F2{T,,Z''^{n),T,0] 



E, 



Ff(^'-)(rfc,z*'^(rfc),r,c) 



by applying definitions (j6.3p and (|6.6p . Thus, 



E 



t,z 



l|„=Tae-'-(^^-*)F2(t,z,r,C)|^T, 



E, 



t,2 



E, 



t,z 



t,z 



E, 



l{,=Tae-^^'^'^-*^i^2(t,^,r,e)|^T, 



l{,=T,}e-''(^'^-*)Ej;^[F2(t,^,r,e)|J^Tj 



F^'(^'-)(Tfc,z(rfc),r,o 



By the preceding identity, (j6.7p and ()6.12p . the identity (|6.9p yields 

M*(t, Z 



sup |e^'' [l|.;,<r^}Fi(t,z,Tfc,ry) 

„<Tfc}+51{„=Tfe} 

F^?(^'-)(rfc,Z(rfc),T,0]]} 



sup <^ E^^ [l|^<2.^}Fi(t, Tfc, r/) 



+1 



sup E, 



Tfc,Z(Tfc) 



F5(^'-)(rfc,z(rfc),r,e)' 



Using the definition (|1.37p of n*, we have 



u*{Tk,z{Tk))= sup Ej'"'^(^'=UF3(^'-)(rfe,z(rfe),r,e) 



and so it follows that 



u*{t,z)= sup l{^^T,,}Fi{t,z,Tk,ri) + lt^r,=Tk}e 



-'•(Tfe-i),,* 



u*{Tk,Z{Tk)) 



Notice that, by the definitions (|6.2p of F"^ and (|6.5p of Fi, we have 

F"*(^'-)(t,z,rfc,r?) = l|,<r,}Fi(t,z,rfc,r?) + ^(Tfc)). 
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The preceding two identities yield 



u*{t,z)= sup E^^ [F"*(^'-)(t,z,rfe,77) 
= u{t,z), (bydSaD). 

This conchides the proof of the theorem. □ 
Proof of Theorem \1.28\ We omit the proof as it is very similar to the proofs of Theorems 11.271 and 

[rT3i □ 

Appendix A. Local a priori boundary estimates 

To complete the proof of Theorem 11.231 we need the following local a priori boundary estimate 
(a parabolic analogue of Corollary 6.7]) for a solution to a parabolic terminal/boundary value 
problem and for which we were not able to find a suitable reference in the literature. 

Proposition A.l (Local a priori boundary estimates). Let ^ M be a domain such that the 
boundary portion Fi is of class C^+". For zq G Fi and R > 0, let 

Br{zo) := i^z eM.'^ : \z - zo\ < and Qr^t{zo) ■= {0,T) x n Br{zo). 

Assume B2r{zq) d H and let f G C°'{Q2r,t{zo)) and g G C'^~^'^{Q2R,t{zq)). Then, there is a 
positive constant C, depending only on zq, R and the coefficients of A, such that for any solution 

U G C2+"(Q2fl,T(2o)) to 

-ut + Au = f on Q2R,t{zo), 

u = g on ((0, T] x {B2r{zo) D Fi)) U ({T} x ^ n B2r{zo)) , 

we have 

II'"IIc2+"{Qk,t(^0)) ^ ^ (ll/lb"(Q2ii,T{^0)) + \\9\\c^+'-{Q2R,Tizo)) + \Mc(Q2R,Tizo))) ' 

Proof. The result follows by combining the global Schauder estimate \34:\ Theorem 10.4.1] and the 
localization procedure of \34:\ Theorem 8.11.1], exactly as in the proof of |16| Theorem 3.8]. □ 

Remark A. 2. The interior version of Proposition lAT] can be found in |34| Exercise 10.4.2]. 

Appendix B. Regular points and continuity properties of stochastic 

representations 

For the purpose of this section, we let d he a non-negative integer, D C M.'^ a bounded domain 
and ti < t2- We denote by Q := (ti,t2) x D and recah that dQ := (ii,i2) x dD U {t2} x D. We 
consider coefficients o, b and a satisfying the following conditions. 

Hypothesis B.l. Let 

a:Q^ M''^'^ and b : Q ^ R'^, 

be maps with component functions, a*-', 6*, belonging to C^'^{Q). Require that the matrix, a, be 
symmetric and obey 

d 

Y,a'^{t,z)ee >m', VCGM^ V(t,z)GQ, (B.l) 
«,i=i 

where (5 is a positive constant. □ 
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Let o" be a square root of the matrix a such that a G C^'^[Q;]R.'^^'^). Such a choice exists by \24\ 
Lemma 6.1.1]. We consider an extension of the coefficients b and a from Q to M x M*^, such that 
these extensions are bounded and uniformly Lipschitz continuous, and condition (jB.ip is satisfied 
on M X M.^. Then, by [30l Theorems 5.2.5 & 5.2.9], for any {t, z) G M x M'^, there is a unique strong 
solution to 

d 

dZi{s) = b'{s,Z{s))dt + '^a'^{s,Z{s))dWj{s), Vi = l,...,d, s > t, 

j=i (B-2) 

Z{t) = z, 

where is a M'^-valued Brownian motion. 

We have the following consequence ofl \12\ Theorem 2.4.2] and \12\ Theorem 2.4.1 and the 
Remark following Theorem 2.4.1]. 

Corollary B.2 (Continuity of stochastic representations with killing term). Assume Hypothesis 
\B.1\ holds and 

(1) the function g is a Borel measurable, bounded function on dQ which is continuous at {t, z), 

(2) the function c : Q ~^ [0; oo) is non-negative, bounded and Borel measurable, 

(3) if there is T > 0, such that tq < T a.s., then the function c : Q ^ M. is bounded and Borel 
measurable function. 

Then 



lim E 

Q3{t' ,z')^(t,z) ^ 

for all regular points [t, z) G <5Q. 



exp ( - / c(s, Z{s))ds ) g{TQ, Z{tq)) 



g{t,z), (B.3) 



Proof. We consider first the case when the stopping time tq is not necessarily bounded by a positive 
constant T. Then, we let cq be a positive constant such that 

< c < Co, a.e. on Q. (B.4) 

Let {t, z) G (5Q be a fixed regular point. We fix e > and consider t' G [ti, ^2] such that \t—t'\ < e/2. 

t' z' 

Then, using the fact that Tq > t' > t — e/2, we see that 



Tq'""' < i - e| C |t - e/2 < t' < r^'"^' < t - e| = 0, 



and so, we obtain 

'1 t',z' 



t\>e]c {rj"' >t + e]u {r'^'^'t < t - e] 



Applying [12i Theorem 2.4.1 and the Remark following Theorem 2.4.1], with tQ := t + e, it follows 
lim Q*''^' {\tq -t\>e)< lim Q*''-'' {tq > t + e) = 0, 

QBit',z')^it,z) VI V I Q9(t/^^/)^(t^2) 

t' z' 

from where it follows that Tq' converges in probability to 0. Similarly, we can argue that 



_t',z' 



exp - / c{s,Z^'''''\s))ds (B.5) 



''See also [SO] Theorem 4.2.12]. 
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converges in probability to 1, as {t',z') G Q tends to {t,z). We again fix e G (0,1) and consider t' 
such that |t' — t\ < — l/(2co) log(l — e). By inequality ()B.4p . we see that 



< 



exp 



> £ 



c{s, Z{s))ds 

exp ^— c(s, Z{s))ds^ < 1 — £^ , (as c > 0), 

[exp {-co{tq - t')) < 1 - e) , (as < c < cq), 

TQ>t' - — log(l - e) 
Co 



TQ>t - ^ log(l - e) ) (because \t' - t\ < -l/(2co) log(l - e)). 
ZCq J 

Choosing to '■= t — log(l — e)/(2co) in \12\ Theorem 2.4.1 and the Remark following Theorem 
2.4.1], we see that the last term in the preceding sequence of inequalities converges to 0, and so 
the collection of random variables (IB.SP converges in probability to 1, as {t' , z') € Q tends to (t, z). 
The sequence is uniformly bounded by the constant 1, and so |20l Exercise 2.4.34 (b)] implies that 
the sequence converges to 1 in expectation also, that is 



lim E: 
From the sequence of inequalities. 



t'z' 



exp 



c(s, Z{s))ds 



1 



0. 



(B.6) 



E, 



t'z' 



exp 



c{s,Z{s))ds g{TQ,Z{TQ)) 



< 



< 



E, 



t',z' 



[9(.tq,Z{tq))] -g{t,z] 



E, 



t',z' 



git,z 

1 — exp 



JEq"" [g{TQ,Z{TQ))]- g{t,z) + ||5||icx,(gQ)E, 



t',z' 



1 — exp 



c{s,Z{s))dsj j g{TQ,Z{TQ)) 
c{s, Z{s))ds 



the conclusion (1B.3P follows from (1B.6P and |12l Theorem 2.4.2] which shows that 

^ u .^q'"' i9{TQ,Z{TQ))]= g{t,z). 

QB(t' ,z')^{t,z) ^ 

We next consider the case when the stopping time tq is bounded a.s. by a positive constant T. 
We fix (t, z) G dQ. Without loss of generality, we may assume that t G [0, T] and Q C [0, T] x M*^. 
Because c is a bounded function on Q, we let ci, C2 be two positive constants such that 



and we set 



and 



-ci < c < C2 a.e. on Q, 
c := c + ci on Q, 



~g{t',z'):=e'^^''-'^g{t',z'), y{t',z')edQ. 

Notice that c is a non-negative, bounded Borel measurable function on Q. Also, 5 is a bounded, 
Borel measurable function on SQ, and it is continuous at (t, z) with 

~g{t,z)=g{t,z). (B.7) 
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In addition, we have for all (t', z') G Q, 

exp {- j^J c{s, Z''^^'{s))ds^ g{TQ, Z'''^'{tq)) 



= exp (^-j^^'' c{s,Z''^^'{s))ds^ ~g{TQ,Z''^^'{TQ)) (B.8) 

+ (exp (ci(t - t')) - 1) exp (- l^^"" cis, ^(tq, ^'(tq)). 

The functions c : Q — )• [0, oo] and g : (5Q — )• M satisfy the requirements of the preceding case, and 
so, we have that 



lim e!1'^' 



exp ( - / c(s, Z{s))ds ] g{TQ, Z{tq)) 



Q5{t' ,z')-^{t,z) ^ 

using (|B.7p . By the boundedness of c on Q, of g on SQ, and the fact that tq <T a.s., we also have 



hm Ei!"^' 



QB{t' ,z')^{t,z) 



(exp (ci(t - t')) - l) exp ^- c(s, Z(s))ds^ g{TQ, Z{tq)) 



0. 



Therefore, the conclusion of the corollary follows from the preceding two limits and identity (jB.Sp . 

□ 
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